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We study the generation of large scale primordial magnetic fields by a cosmological phase tran- 
sition during the radiation dominated era. The setting is a theory of N charged scalar fields cou- 
pled to an abelian gauge field, that undergoes a phase transition at a critical temperature much 
larger than the electroweak scale. The dynamics after the transition features two distinct stages: 
a spinodal regime dominated by linear long-wavelength instabilities, and a scaling stage in which 
the non-linearities and backreaction of the scalar fields are dominant. This second stage describes 
the growth of horizon sized domains. We implement a formulation based on the non-equilibrium 
Schwinger-Dyson equations to obtain the spectrum of magnetic fields that includes the dissipative 
effects of the plasma. We find that large scale magnetogenesis is efficient during the scaling regime. 
Charged scalar field fluctuations with wavelengths of the order of the Hubble radius induce large 
scale magnetogenesis via loop effects. The leading processes are: pair production, pair annihilation 
and low energy bremsstrahlung, these processes while forbidden in equilibrium are allowed strongly 
out of equilibrium. The ratio between the energy density on scales larger than L and that in the 
background radiation r(L,T) = ps(L, T)/p crn b(T) is r(L,T) ~ 10 -34 at the Electroweak scale and 
r(L,T) ~ 10~ 14 at the QCD scale for L ~ 1 Mpc. The resulting spectrum is insensitive to the 
magnetic diffusion length and equipartition between electric and magnetic fields does not hold. We 
conjecture that a similar mechanism could be operative after the QCD chiral phase transition. 
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I. INTRODUCTION 

A variety of astrophysical observations including Zeeman splitting, synchrotron emission, Faraday rotation mea- 
surements (RM) combined with pulsar dispersion measurements (DM) and polarization measurements suggest the 
presence of large scale magnetic fields [l|, |2|, y, |4j, yj, M ■ The strength of typical galactic magnetic fields is measured 
to be ~ \x G0Ty, 0, @ and they are correlated on very large scales up to galactic or even larger reaching to scales 
of cluster of galaxies ~ 1 Mpcg, H H @. The origin of these large scale magnetic fields is still a subject of much 
discussion and controversy. It is currently agreed that a variety of dynamo mechanisms are efficient in amplifying 
seed magnetic fields with typical growth rates T ~ Gyr -1 over time scales ~ fO — 12 Gyr (for a thorough discussion 
of the mechanisms and models seeQ). The ratio of the energy density of the seed magnetic fields on scales larger 
than L (today) to that in the cosmic background radiation, r(L) = pb{L)/ p cm b must be r(L ~ IMpc) > 1CP 34 for a 
dynamo mechanism to amplify it to the observed value, or r(L ~ IMpc) > 10~ 8 for the seed to be amplified solely 
by the gravitational collapse of a protogalaxypj. 

There are also different proposals to explain the origin of the initial seed. Astrophysical batteries rely on gradients 
of the charge density concentration and pressure and their efficiency in producing seeds of the necessary amplitude 
is still very much discussed [3, |4|. Primordial magnetic fields that could be the seeds for dynamo amplification can 
be generated at different stages in the history of the early Universe, in particular during inflation, preheating and 
or phase transitions 0, 0, 0. Primordial (hyper) magnetic fields may have important consequences in electroweak 
baryogenesis|7|, Big Bang nucleosynthesis (seej3j), the polarization of the CMB[8| via the same physical processes 
as Faraday rotation, and structure formation^ 0, || , thus sparking an intense program to study the origin and 
consequences of the generation of magnetic fields in the early I Jniverse [lfij - pfij . 

A reliable estimate of the amplitude and correlations of seed magnetic fields must include the dissipative properties 
of the plasma, in particular the conduct ivitv[Tll fiH, fl4| . In ref.[27| we have introduced a formulation that allows to 
compute the generation of magnetic fields from processes strongly out of equilibrium. This formulation, which is based 
on the exact set of Schwinger-Dyson equations for the transverse photon propagator is manifestly gauge invariant and 
is general for any matter fields and any cosmological background (conformally related to Minkowski space-time). In 
the case in which strongly out of equilibrium effects arise from long-wavelength fluctuations, such as during phase 
transitions, this formulation allows to separate the contribution of the hard degrees of freedom which are in local 
thermodynamic equilibrium from that of the soft degrees of freedom that fall out of LTE (local thermal equilibrium) 
during the phase transition and whose dynamics is strongly out of equilibrium. This separation of degrees of freedom 
leads to a consistent incorporation of the dissipative effects via the conductivity (for details seejllj). In that reference 
a study of magnetogenesis in Minkowski space-time during a supercooled phase transitions was presented and the 
results highlighted the main aspects of the generation of magnetic and electric fields in these situations. 

We study the generation of large scale (hyper) magnetic fields by a cosmological phase transition during a radiation 
dominated era by implementing the formulation introduced in ref.|27|. The setting is a theory of N charged scalar 
fields coupled to an abelian gauge field (hypercharge). We consider the situation when this theory undergoes a phase 
transition after the reheating stage and before either the Electroweak or the QCD phase transition, since we expect 
that these transitions will lead to new physical phenomena. The non-perturbative dynamics out of equilibrium is 
studied in the limit of a large number TV of (hyper) charged fields and to leading order in the gauge coupling. The 
non-equilibrium dynamics of the charged scalar sector features two distinct stages. The first one describes the early 
and intermediate time regime and is dominated by the spinodal instabilities which are the hallmark of the process of 
phase separation and domain formation and growth. This stage describes the dynamics between the time at which 
the phase transition takes place and that at which non-linearities become important via the backreaction. The second 
stage corresponds to a scaling regime which describes the slower non-equilibrium evolution of Goldstone bosons and 
the process of phase ordering|26j and growth of horizon-sized domains. This scaling regime is akin to the solution 
found in the classical evolution of scalar field models with broken continuous symmetries after the phase transition 
that form the basis for models of structure formation based on topological defect sp8H29j. 

The solution of the scalar field dynamics [2(| is the input in the expression for the spectrum of the magnetic field 
obtained in |27j| to obtain the amplitude of the primordial seed generated during both stages. 

We find that scaling stage is the most important for the generation of large scale magnetic fields. Large scale 
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magnetic fields are generated via loop effects from the dynamics of modes that are at the scale of the horizon or smaller. 
The leading order processes that result in the generation of large scale magnetic fields are: i) pair production, ii) 
pair annihilation and iii) low energy bremsstrahlung. These processes would be forbidden in equilibrium by energy 
momentum conservation, but they are allowed strongly out of equilibrium because of the rapid time evolution of the 
cosmological background and the fast dynamics of the scalar field fluctuations. 

The resulting spectrum is rather insensitive to the diffusion length scale which is much smaller than the horizon 
during the radiation dominated era. The ratio of the magnetic energy density on scales larger than L (today) to the 
energy density in the background radiation r(L, 77) = pb{L,t])/ p cm b(v) is summarized in a compact formula [ea. (|5.99|) ]. 
For L ~ 1 Mpc (today) we find r(L,ri) ~ 10 -34 at the Electroweak scale and r(L,ri) ~ 10~ 14 at the QCD scale, 
suggesting the possibility that these primordial seeds could be amplified by dynamo mechanisms to the values of the 
magnetic fields consistent with the observed ones on these scales. 

II. THE PHYSICAL PICTURE 

Whether a symmetry breaking phase transition occurs in local thermodynamic equilibrium (LTE) or not is a 
detailed question of time scales. Two important time scales must be compared, the time scale of cooling through the 
critical temperature t coo ; = T(i])/T(ri) and the relaxation time scale of a fluctuation of typical wavelength k, r(k). 
If cooling through T c occurs on time scales slower than the relaxation scale, i.e, r coo / 3> T(k) then a fluctuation of 
wavevector k adjusts to LTE and for these scales the phase transition occurs in LTE. On the other hand if r coo ; <C r(fc) 
fluctuations of wavelength fc" 1 cannot adjust to the conditions of LTE and fall out of equilibrium, i.e, freeze out and 
for these scales the phase transition occurs on very short time scales like a quench from a high temperature phase 
into a low temperature phase [3(| • While short wavelength modes typically have very short relaxation time scales and 
remain in LTE during the phase transition, long-wavelength modes undergo critical slowing down near T c and their 
relaxational dynamics becomes very slow[3l|. Thus as the temperature nears the critical with a non- vanishing cooling 
rate, long-wavelength modes freeze out falling out of LTE, namely for these modes the phase transition occurs out of 
equilibrium. In theories in which the equilibrium phase transition is of second order, the non- equilibrium dynamics 
below the critical temperature is described by the process of spinodal decomposition: long-wavelength fluctuations 
become unstable and grow. In Minkowski space-time this growth is exponential in time but in an expanding cosmology 
the growth will depend on the time dependence of the scale factor. 

The growth of long-wavelength fluctuations entails that the field becomes correlated within regions characterized 
by a time-dependent correlation length £(7/) and the amplitude of the long-wavelength fluctuations becomes non- 
perturbatively large. The non-perturbative amplitudes of long-wavelength fluctuations, becoming of 0(1/ ^/\) with 
A the scalar self-coupling are a consequence of the fact that the mean-square root fluctuation of the field probes the 
broken symmetry states with (</> 2 ) oc 1/A with A. 

This is the hallmark of the process of phase separation, the correlated regions correspond to domains, inside 
these domains the field is near one of the vacuum states. A particular time scale, the nonlinear time t n i tx ln(l/A) , 
determines the transition from a regime of linear instabilities to one in which the full non-linearities become important. 
This time scale roughly corresponds to when the mean square root fluctuation of the scalar field samples the broken 
symmetry states and the phase transition is almost complete. At this stage the amplitude of the long-wavelength 
fluctuations become of C(l/A), the phases freeze out and the long- wavelength fluctuations become classical but 
stochastic^- (33J. 

Consider now the situation where the scalar field carries an abelian charge and is coupled to the electromagnetic 
gauge field. The strong spinodal fluctuations in the scalar field will induce fluctuations in the current-current correla- 
tion function, and while the expectation value of the current must vanish by translational and rotational invariance, 
the current and charge correlators will have strong non-equilibrium fluctuations. These current fluctuations will in 
turn generate a magnetic field with a typical wavelength corresponding to the wavelength of the spinodally unstable 
modes. 

This is the main premise of this article: the spinodal instabilities which are the hallmark of a non-equilibrium 
symmetry breaking phase transition will lead to strong charge and current fluctuations of the charged scalar fields 
which in turn, lead to the generation of magnetic fields through the non-equilibrium evolution. 

The main ingredients that must be developed in order to understand the generation of magnetic fields through this 
non-equilibrium process are: 

• A consistent framework to compute the spectrum of generated magnetic field, namely (B(k, t) ■ B(—k, t)} /V with 
B(k,t) the spatial Fourier transform of the Hcisenbcrg magnetic field operator and V the (comoving) volume of 
the system. 
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• We anticipate that plasma effects must be included to assess the generation and eventual decay of the magnetic 
fields. If there is a large conductivity in the medium the magnetic field will diffuse but also its generation will 
be hindered. This will be a point of particular importance within the cosmological setting pi IfiL fill Il3l Il4| . 

• A major challenge of any mechanism of large scale magnetogenesis is to generate the seed magnetic fields 
from microscopic, causal processes. An important aspect of the results presented here is that this generation 
mechanism is mediated by loop effects and correspond to processes that are forbidden in equilibrium but allowed 
strongly out of equilibrium. 

These ingredients will be analyzed in detail below. 



III. MAGNETIC FIELDS IN FRIEDMANN-ROBERTSON- WALKER COSMOLOGY 



The cosmological setting in which we are primarily interested corresponds to a symmetry breaking phase transition 
in a radiation dominated Universe. Such phase transition is in principle different from the electroweak one 1 and 
presumably occurs at a much higher energy scale, such as the GUT scale ~ 10 15 Gev but is assumed to be described 
by a particle physics model that includes many fields with (hyper)-charge either fcrmionic or bosonic. We will not 
attempt to study a particular gauge theory phcnomenologically motivated by some GUT scenario, but will focus our 
study on a generic scalar field model in which the scalar fields carry an abelian (hyper)charge. The simplest realization 
of such model is scalar electrodynamics with N charged scalar fields <p r , r = 1, . . . , N and one neutral scalar field ip 
whose expectation value is the order parameter associated with the phase transition. The neutral field is not coupled 
to the gauge field and its acquiring an expectation value does not break the U(l) gauged symmetry. This guarantees 
that the abelian gauge symmetry identified with either hypercharger or electromagnetism is not spontaneously broken 
to describe the correct low energy sector with unbroken U(\)em- We will take the neutral and the N complex 
(charged) fields to form a scalar multiplct under an 0(2N + 1) isospin symmetry. The electromagnetic coupling 
explicitly breaks the 0(2N + 1) symmetry down to SU(N) x U(l). In the absence of electromagnetic coupling as 
the neutral field acquires an expectation value the isospin symmetry is spontaneously broken to 0(2N). Since by 
construction only the neutral field acquires a non vanishing expectation value under the isospin symmetry breaking 
the photon remains massless (it will obtain a Debye screening mass from medium effects, see discussion below). 

The action that describes this theory in a general cosmological background is given by 




AN \ 2 



-.J-ilv -Tap 9 9 



(3.1) 



and 



i_i ■ O^i iGvA/j, and J~ ^ v - d^^Ai/ df-A^ Z6[ w 4^ fi , ^4.^] 



(3.2) 



A' 



N 



Furthermore, anticipating a non-perturbative treatment of the non-equilibrium dynamics of the scalar sector in a 
large TV expansion, we have rescaled the quartic coupling in such a way as to display the contributions in terms of 
powers of 1/N, keeping A fixed in the large N limit. 
A Friedmann-Robertson- Walker line element 

ds 2 = dt 2 -a 2 {rj) dx 2 , (3.3) 
is conformally related to a Minkowski line element by introducing the conformal time r\ and scale factor C(rj) as 

dt 



V 



a(rj) 



; C{rj) = 0(77(77)) 



(3.4) 



If the electroweak phase transition is weakly first order, nucleation will be almost indistinguishable from spinodal decomposition and 
the phenomena studied here may be of relevance. 
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In terms of these the line element and metric are given by 

ds 2 = C 2 (n) (dn 2 -dx 2 ) ; 9llv = C 2 (n) Vfw , (3.5) 
where rj^u = diag(l, — 1, — 1, —1) is the Minkowski metric. Introducing the conformal fields 

Apfji, x) = A^tirj), x) ■ $(77, x) = C(r?) (f>(t(r,),x) ; x) = C(r?) ip(t(r)), x) 
and in terms of the conformal time, the action now reads 



S= I drj d 3 x 
with 



(3.6) 



C"(n) 

M 2 (r ] )=- f j, 2 C 2 (r ] )--^- , Dp = dp-ieA ll ; = d»A v - d v A„ , (3.7) 



and the primes refer to derivatives with respect to conformal time. Obviously the conformal rescaling of the metric 
and fields turned the action into that of a charged scalar field interacting with a gauge field in flat Minkowski space- 
time, but the scalar field acquires a time dependent mass term 2 . In particular, in the absence of electromagnetic 
coupling, the equations of motion for the gauge field A^ are those of a free field in flat space time. This is the 
statement that gauge fields are conformally coupled to gravity and no generation of electromagnetic fields can occur 
from gravitational expansion alone without coupling to other fields or breaking the conformal invariance of the gauge 
sector. The generation of electromagnetic fields must arise from a coupling to other fields that are not conformally 
coupled to gravity, or by adding extra terms in the Lagrangian that would break the conformal invariance of the gauge 
fields [13 ■ 

The conformal electromagnetic fields £, B are related to the physical E, B fields by the conformal rescaling 

£ ft 

E= c 2 ~W) 5 B =c 2 ~W)' (3 - 8) 

corresponding to fields of scaling dimension two. 

Gauge invariance: Since we will study the generation of magnetic fields through loop effects and keeping only low 
order terms in the diagrammatic expansion, it is important to guarantee that the results obtained are gauge invariant, 
namely physical. 

A gauge invariant formulation leads to the following Lagrangian density (for details see pTjl^ ) 

C = l -3^ + 0^ + l -d^A T ■ d»A T + i(VAo) 2 - M 2 ( v ) Q* 2 + $t$^ _ Q* 2 + $t<j, x 

-iel T • ($ f V$ - V$ f $) - e 2 (A| - A%) $ f <I> - ie A a - , (3.9) 

where $ is a gauge invariant local field which is non- locally related to the original fields, and At is the transverse 
component of the vector field (V • At = 0) and Ao is a non-propa gati ng; field as befits a Lagrange multiplier, its 
dynamics is completely determined by that of the charge density [271 l34j . 

The main point of this discussion is that the framework to obtain the power spectrum of the generated magnetic 
field presented below is fully gauge invariant. 

IV. PHASE TRANSITIONS IN RADIATION DOMINATED COSMOLOGY 

A. Kinematics 

We consider here a phase transition in a radiation dominated cosmology in a theory with a critical temperature 
near but below a GUT scale, and we take the reheating temperature of the post-inflationary universe to be higher 



2 Here we neglect the effect of the conformal anomaly |20l 



6 



than the critical temperature. Furthermore we assume that after reheating the Universe is described by a state of 
local thermodynamic equilibrium and the background radiation dominates the dynamics of the Hubble expansion. 

In this radiation dominated cosmology, as the initial state local thermodynamic equilibrium at an initial temperature 
Tr > > T c the expansion results in cooling and the temperature eventually falls below the critical triggering the phase 
transition. 

Using finite temperature field theory in an expanding background geometry, it is shown pfij that the effective time 
dependent mass term depends on the effective time dependent temperature T(rj) = Tji/a{rj) (see discussion below) 
which reflects the cooling from the cosmological expansion. Hence at a given time the temperature equals the critical 
and the phase transition occurs. Field modes with wavectors much larger than the symmetry breaking scale \x will 
remain in LTE and will not be affected by the symmetry breaking dvnamics|30|. but those with wavevectors much 
smaller than the symmetry breaking scale will fall our of equilibrium and become spinodally unstable |30|. 

We normalize the scaling factor C{rf) at the reheating time rj — t\r in such a way that C(j]r) = 1 then the explicit 
expression for C(rj) reads 

C( V )=H R r, (4.1) 

where Hr is the Hubble constant at the reheating time, Hr = r]^ 1 . 

We can relate Hr to the reheating temperature and the Planck mass G~ 5 through the equation 



30 

and the Einstein-Friedman equation 



n 9 * T R (4.2) 



1/2 T 2 



where g* is the effective number of degrees of freedom at the reheating temperature and we introduced the scale M* 
of the order of the Planck mass 

M» = 3 -4 * . (4.4) 
2 7T2 V3* " 

In radiation dominated epoch the time-dependence of the mass term (I3.7|l is given by the expression 

-M\n) = H 2 R ?7 2 = mV, (4.5) 
where we see the emergence of a new mass scale 



/~i = y/ji Hr . (4.6) 

This scale will play an important role in the following discussion and in the comparison with results obtained in 
Minkowski space-time [27|]. There is a last scale which plays a relevant role, the horizon scale rn(jl) which is fixed by 
the evolution on the time of the Hubble constant: 

r H (v) = j^tt = C( V ) n = H R r? . (4.7) 
Modes with physical wavelength X P h vs = i, 2lT inside the horizon 

" y fcphys 

X phys (v) - k~h ys (V) < r H fa) (4.8) 

are causally connected, modes outside the horizon are causally disconnected. 
The relaxation rate of hard modes of the charged fields is given bv|35| 

r(n) - a T{n) In - , (4.9) 
a 



where the effective temperature varies with time as 



T(V) - (4.10) 
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and the expansion rate given by 



T 2 (ri) 

H(n) = — . (4.11) 



Therefore, 

r(r?) 10 16 



H( V ) T(r?)[Gev] 



(4.12) 



Thus hard modes are in thermal equilibrium for T(ry) < 10 15 Gev. 

In particular, modes with k ~ Tr are the hard modes that give the leading contribution to the conductivity in 
the high temperture limit [3ft l37j. Modes with k < u will manifest the long- wavelength spinodal instabilities and 
their dynamics will be strongly out of eauilibrium|30L l32l Li'l |3S| . Their amplitude becomes non-perturbatively 
large j30t 1321 1331 138| and will be responsible for the non-equilibrium generation of the primordial magnetic field j27|. 

Using eqs. ||0]| and l|4.1U|l we can write the conformal time as 

Tr M. 

^=H R ^ = T R -m- (4 - 13) 

As it will become clear below an important cosmological quantity is the product 

u k r>t \ k phys(i) rr, s / \ 27r M * u ,„\ 

k 71 = cW) C{71) 71 = ~rW iv) rHiv) = l7F r tW)- ( } 



The ratio 

kphys(v) _ 27T 



(4.15) 



is a kinematical invariant. Its value today is determined by the scale L which will be typically chosen to correspond 
to a galactic scale or the scale of galaxy clusters, and the temperature of the CMB. It is given by, 



lt ^ 3 - 7 x 1025 (m^) 



(4.16) 



Therefore, 



1W /Mdc\ ( 10_22 for T ^ = T R - 1015 Gev 
k r] ~ 10" 9 -f^ -f^ = I 10" 9 for the EW transition (4.17) 

T ^> V L J y io- 6 for the QCD transition 

for L ~ 1 Mpc. Thus, during the regime of interest in this article, kr\ <C 1 for scales of galaxy clusters. A noteworthy 
aspect of ea. l|4.17|) is that the wavelengths corresponding to the scale of galaxies or clusters today were well outside 
the horizon during the radiation dominated era when the electroweak and QCD phase transitions occurred. 

Thus the challenge for any mechanism of primordial magnetogenesis operating on scales well inside the Hubble 
radius, such as causal microscopic processes, must generate magnetic fields with wavelengths which are much larger 
than the Hubble radius. 

Another important quantity is the ratio of the wavevector k of the primordial magnetic field to the conductivity. 
As it will be discussed below, the physical conductivity is given by 

C N( V ) T( V ) 

<r{V) = ] 1 i (4-18) 

where C is a constant of 0(1), N(rj) is the number of ultrarelativistic charged species, and we have neglected the 
(logarithmic) dependence on the energy scale in the running coupling constant. For this discussion we will neglect 
the time dependence of N(rj) assuming that the number of charged ultrarelativistic species remains constant (this 
assumption can be relaxed without qualitative modifications of the main argument). Under this assumption 



(4.19) 
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with or being the comoving conductivity determined at the time of reheating 



(4.20) 



Thus the ratio, 



neglecting logarithmic corrections. 
Furthermore, 



kphysiv) 27ra _ 10 _ 27 /Mpc\ 



o(rj) NLT R \ L 



10 5 for T(r?) = Tr ~ 10 15 Gev 



d{ri)r H {ri) ^ "( 1018 for the EW phase transition (4.22) 

aT W) I 10 2i for the QCD phase transition 

where we have neglected logarithmic corrections. Therefore or i] 3> 1 throughout the radiation dominated era 
considered in this article. The regime or r\ 3> 1 ; k 2 rj/ 'or <C 1 is dominated by the (slow) hydrodynamic relaxation 
of the magnetic field. 



Another relevant estimate involves the (comoving) diffusion length £,<uff(jl) = \Zv/ (J R 




for T{rj) =T R ~ 10 15 Gev 

for the EW phase transition (4.23) 
for the QCD phase transition 

where again we have neglected logarithmic terms. Therefore the diffusion length is much smaller than the Hubble 
radius during the radiation dominated era. Finally, combining (|4.23|) and (|4.17|l we find 

1(T 25 < kUffiv) < lO" 16 , (4.24) 

between reheating and the time of the QCD phase transition. 

The contribution from the hard modes of both the charged scalar and gauge fields which remain in local ther- 
modynamic equilibrium lead to an effective mass for the scalar field. This thermal mass is obtained from the long- 
wavelength limit of the scalar field self-energy and includes the hard thermal loop contributions from the gauge and 
scalar fields |42ll43|. This thermal mass is given by 



m 



»=||(A + 3,;) . (-1.25) 



Finally, another important quantity is the Debye screening length that determines the scale at which long-range forces 
are screened by the polarizability of the medium. In an ultrarelativistic plasma, the comoving Debye screening length 
is given by|H Hj 

£d ~ 4fr (4-26) 
e T R 

the ratio of the Debye screening length to the Hubble radius is given by 

Z-H^lltnl (4 27) 

Hence £d "C for T(rf) < 10 16 Gev, thus long range forces are screened over very short distances. The formation 

of long-wavelength domains with typical size of the order of the Hubble radius 26] leading to strong charge and current 
fluctuations that will seed magnetic fields, will not be hindered by long-range forces, which are effectively screened 
over sub-horizon distances. 

Magnetic field generation via charge asymmetries during a period in which electromagnetism was spontaneously 
broken was previously studied by Dolgov and Silk|3 who argued that long-range forces would be screened by the 
Higgs mechanism. This is different from the situation studied in this article, where the U (1) symmetry associated with 
electromagnetism (rather hyper charge) is not spontaneously broken. Long range forces are screened by the plasma, a 
situation not considered inl44l . 



9 



B. Scalar fields dynamics 



For completeness and to highlight the aspects of the non-equilibrium dynamics most relevant to the generation of 
magnetic fields, we summarize the main features of scalar field dynamics. For further details the reader is referred 
to [32 . l33l l38| . In what follows we will neglect the backreaction of the gauge fields on the dynamics of the scalar fields. 
The rationale for this is that the main non-equilibrium processes that lead to magnetogenesis will be non-perturbative 
in the scalar sector and result from the instabilities associated with the phase transition. The contribution from 
the gauge fields, in the form of self-energies for the scalar fields, do not feature the instabilities associated with the 
phase transition and will, furthermore, be suppressed at least by one power of a the (hyper) electromagnetic coupling- 
constant as compared to the scalar self-interaction. 

As described above, the non-equilibrium evolution of long-wavelength modes begins with the spinodal instabilities 
which result in an exponential growth of the amplitudes for long-wavelength fluctuations. When the non-linearity 
becomes of the same order as the tree-level terms in the equations of motion, the back reaction of these fluctuations 
shuts-off the instabilities |32l 1.3.1 l.3q . Therefore a non-perturbative treatment of the dynamics is required. The large 
TV limit of the scalar sector allows a systematic non-perturbative treatment of the dynamics which is rcnormalizable 
and maintains the conservation laws|32ll38|. 

We will therefore study the dynamics in leading order in the large N limit that already reveals the important 
non-equilibrium features of the evolution. 

The contribution from the gauge fields to the equations of motion of the long- wavelength modes of the scalar fields 
arise through self-energy corrections. To lowest order in a these are dominated by the hard modes of the gauge fields 
(hard thermal loops) which lead to a contribution to the thermal mass given by eT/V8 jHEJ and had already been 
accounted for in the thermal mass (|4.25() . 

The non- equilibrium effects in the gauge contribution of the scalar self-energy will arise from polarization loops in 
the intermediate photon lines and are, therefore, of 0(a 2 ). Since our calculation will be to lowest order in a, we will 
neglect these contributions. Thus, to this order the gauge field contribution to the scalar self-energy is accounted for 
in the thermal mass. 

Hence the dynamics of the scalar field is studied along the same lines as presented in refs. [32LI33 . 38] but the only 
difference is in the initial conditions in the modes that reflect the thermal mass in LTE. 
Since symmetry breaking is chosen along the direction of the neutral field W we write 

(4.28) 



where the expectation value is taken in the time evolved density matrix or initial state. The leading order in the large 
N limit is obtained either by introducing an auxiliary field and establishing the saddle point or equivalently by the 
factorizations [32, |38j 



The non-linear terms of the a field lead to contributions of 0(1/ N) in the large N limit, and to leading order 
the dynamics is completely determined by the N complex scalars $ r . This factorization that leads to the leading 
contribution in the large N limit makes the Lagrangian for the scalar fields quadratic (in the absence of the gauge 
coupling) at the expense of a self-consistent condition: thus charged fields <f> acquire a self-consistent time dependent 
mass. 

The dy namics is determined by the Heisenberg equations of motion of the neutral field and the charged fields 
3> [32, 133L l38j| . We will consider that at the onset of the radiation dominated era, the system is in the symmetric high 
temperature phase in local thermal equilibrium with a vanishing expectation value for the scalar fields. In the absence 
of explicit symmetry breaking perturbations the expectation value of the scalar field will remain zero throughout the 
evolution, thus <p = 0. 

It is convenient to introduce the mode expansion of the charged fields 



d 3 k 



a r (k)h( V ) e^ s + bi(k) f* k (ri) e- 



ik-x 



r = l. 



,N 



(4.29) 



In leading order in the large N limit, the Heisenberg equations of motion for the charged fields translate into the 
following equations of motion for the mode functions for 77 > r]R [3^, IU 



^ + fc 2 -M 2 (,) + ^(,) + A ($ t^ 



f k (v) = . 



(4.30) 
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We must now append initial conditions for the mode functions fk(v)- The initial conditions on the mode functions 
fk(v) depend on the value of the wavevector k as compared to the horizon scale H^ 1 at the reheating time: 

• k > Hr: for fluctuations inside the horizon, we may assume thermal quasi-particle boundary conditions at the 
reheating temperature: 

f k ( m ) = -L= ; f k ( m ) = -iWk f k {m) , w k = Jk 2 + m 2 T , (4 31) 

where the frequencies W k are quasi-particle frequencies with thermal mass m\ given by eq. 14.25( 1 at the reheating 
temperature, 



W k = J k 2 + m 2 T . 

For these modes the assumption of local thermodynamic equilibrium is well motivated and we have 

(4(fc) a s (k)) = <&t(fc) b s (k)) = S rs n k ■ n k = -^ (4.32) 

e t — 1 



• k < Hr: for supcrhorizon fluctuations, which are causally disconnected at the reheating time, we cannot assume 
a thermalized distribution. The correct distribution has to be derived by following the dynamics from the 
inflationary stage, when the fluctuations were well inside the horizon. While a complete discussion of the initial 
conditions is left to a forthcoming article, the case under consideration we will see that the dependence on the 
initial conditions is rather weak and only during the initial stages of the phase transition. For the later stages, 
dominated by the scaling solution described below, the dynamics is universal and does not depend on the initial 
conditions. We will simply assume that both W k and n k have a finite non-zero limit as k — > namely the only 
important quantities for the dynamics of long-wavelength fluctuations are 

lim W k = W , < W < oo, lim n k = n . < n < oo . (4.33) 

fe^O fe^O 

With this choice of the initial state we find the backreaction term to be given by 



$t$> = I / l/fc(?7)|2[1 + 2nk] ■ (4 - 34) 



_x_ 

2N^ 4 J (2?r) 

This expectation value is ultraviolet divergent, it features quadratic and logarithmic divergences in terms of an upper 
momentum cutoff. The quadratic divergence and part of the logarithmic divergence (the one proportional to the 
mass term) are absorbed in a renormalization of the mass term /i 2 — ► [i 2 R and the remainder logarithmic divergence is 
absorbed into a renormalization of the scalar coupling A — ► Xr. While these aspects are not relevant for the discussion 
here, they are mentioned for completeness, the reader is referred to j2(| for details. 

After renormalization the self-consistent field t^t is subtracted twice, and is given by (for details see |2(J and 

references therein) 



x A fl r 2j ,,., |2 i e{ q -K 2 ) 

Xr = ^2 J o 9 dq i \f q (r])\ 2 



q 2q 



3 



-/4 + ^ 7 (* t 'i>> 



(4.35) 



and the mass and coupling are replaced by their renormalized counterparts fi R ; Xr respectively. Here K is an arbitrary 
renormalization scale. In order to avoid cluttering of notation we now drop the subscript R for renormalized quantities, 
in what follows /i; A stand for the renormalized quantities. 

The finite temperature term J(rj) has contributions from short wavelengths for which the mode functions are of 
the form f q (ji) ~ e lqn / \fq and contributions from long wavelengths. The contribution from short wavelengths is the 
same as that in equilibrium in Minkowski space time and determines the hard-thermal loop [42l liflj contribution to 
the self-energy given by|2(j 

Jhtl = H (4.36) 
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where we have used that the short wavelength modes are in thermal equilibrium at the reheating temperature Tr. 
This hard thermal loop contribution has been self-consistently accounted for in the thermal mass of the scalar field 

It is convenient to separate the hard thermal loop component (|4.36[) from eg. 14.3511 and define 



<Z 2 dql |/ ? (r?)| 2 (l 



2n q ) - - 

q 



e T R - 1 J 



Q{q~K 2 



2q 3 



AS(7y) 



(4.37) 



After renormalization and in terms of dimensionless quantities, the non-equilibrium dynamics of the charged scalar 
fields is completely determined by the following equations of motion [2ti 132. 1331 l3q , 

dr/ 2 



M 2 ( V ) 



AS(r?) 



/«(»?) = o ; f g (m) 



f q {m) = -iW q f q { m ) 



with W q given by eq. 1)4. 5J) and the effective, (conformal) time dependent mass is given by 

M 2 {r,) = C 2 {r,) /j 2 
24 fi 2 



T 2 

1 R 



C 2 (v) T 2 



A + 3e 2 ' 



(4.38) 

(4.39) 
(4.40) 



The time dependent mass term M 2 (rj) includes the high temperature corrections and clearly displays the cooling 
associated with the expansion in the form of a time dependent effective temperature T e ff(rj) = Tn/C(r)). The 
phase transition occurs at a time r\ c when T e ff{rj c ) — T Cl thus for rj > rj c the effective time dependent mass term is 
M 2 (rj) — M 2 (r]) = —fi 2 B.\ if = —jl 4 rj 2 as given by equations (|4.5|) - (|4.t)|l . 

The full time evolution of mode functions in a radiation dominated cosmology has been studied analytically and 
numerically in detail in ref . prl l38| . Here we highlight the most important features which are necessary ingredients to 
study magnetogenesis. The reader is referred to[2(j for a more comprehensive discussion. 

The are two main dynamical stages in the evolution: 

• Spinodal stage: this is the stage immediately after the phase transition which is dominated by spinodal decom- 
position and the growth of long- wavelength fluctuations |27| . This stage spans the time scale r\ c < r\ < r) n i where 
the non-linear time scale rj n i is determined by (see below) 



AE(?7„i) 



A 4 



(4.41) 



During this stage the back-reaction, determined by the term A £(f?), can be neglected and the dynamics is 
linear. 

Scaling stage: This is a stage in which the non-linearity encoded by the back-reaction term AS(^) are very 
important and compete with the tree level term in the equations of motion. This stage is described by a scaling 
solution of the equations of motion for the modes with small wavevectors and describes the non-equilibrium 
relaxation of long- wavelength fluctuations [2tiL l28l |29| . 



1. Spinodal stage 

After the phase transition but before the non- linear time scale after which the back-reaction becomes important, 
namely for t) c <C rj < rj n i the time dependent mass term is given by Ai 2 — ~jl A r] 2 , and for weak coupling A < 1 we 
can neglect the back- reaction AS (77). The equations of motion for the mode functions during this stage are given by 



d V 2 



~4 2 

H 77 



SM = 



q < v 2 v 



(4.42) 



We note that for T r - 10 15 GeV 



f 1 Vc = 



/A + 3e 2 M* 
24 IT 



1/2 



10 



(4.43) 



12 



and therefore for n > r\ c we are in the regime /2 77 ^> 1 . It is clear that the mode functions f q {rj) will increase 
exponenti ally in the band of unstable wavevectors q < fPn. Eg. 14.421) can be solved exactly in terms of Hermite 
functions 143 



f q {rj) = b q ") 2 H k j^s [p rj) + a q ^ #_i(4 +1 ) (ifir,) 



(4.44) 



where the constants a q and b q are fixed by the initial conditions Ij4.38|l . For fi r\ 1 we can use the asymptotic 
behavior of the Hermite functions 1451. 



H v {z) (2zY 



l + O 



1 



and we find for the mode functions, 



fqW = a q 



(A v)~ 



l + O 



f? rj 2 



(4.45) 



Since the exponentially damped solution becomes negligible the phases of the mode functions f q {rf) freeze, namely, 
they become constant in time and are slowly varying functions of q for long wavelengths. 

This is very similar to the situation in Minkowski space-time, where the mode functions however increase as e^*, 
i.e. much slower. In any case the soft (q — > 0) modes are the most amplified at the end of the evolution, therefore, 
the quantum fluctuations (|4.37(l are dominated by the lower integration bound q — 0. 

We notice that the freezing of the long-wavelength mode functions will play an important role in the discussion 
about the magnetic field generation, since it assures the independence of the final result from the initial particle 
distribution function, except for subleading corrections. 

The physics of the phase transition is essentially the same as in Minkowski space-time [23, HE H3i smce the 
exponential growth of modes in the spinodally unstable band will make the back reaction term A£(r/) begin to grow and 
eventually cancel the term — fi r] 2 in the equations of motion (for rj >> r/ c the effective time dependent temperature 
vanishes) . 

This will happen at a non-linear time scale defined by |3z 



A T,{r] n i) = fi 4 r% 



(4.46) 



Two important aspects are described by rj n f. i) at this time scale the phase transition is almost complete since 
AS(?7 rl ;) = fi A rfci means that A($^<I>) /2N = fs^n 2 ^-, namely the mean square root fluctuations in the scalar field probe 
the manifold of minima of the potential. 

ii) At n ~ rj n i the mean square root fluctuations of the field are of order M 2 (rj n i)/\ probing the vacuum manifold, 
and the non-linearities become very important. The back reaction AS (rj n i) becomes comparable to M 2 (rj n i) and the 
instabilities shut-off. Thus for r/ c < n < r\ n i the dynamics is described by the linear spinodal instabilities while for 
r\ > rj n i a full non-linear treatment of the evolution is required. As it will be discussed below this later stage is 
described by the emergence of a scaling solution. 

For rj n i > i] ^> the asymptotic form (|4.45|l for the mode functions apply and we find for the the quantum 
fluctuations (|4.37l) which dominated by the lower integration bound q = 0, 



AE(77 ni ) = A (1 + 2 n ) 



fi 2 |«o| 2 



-2 2 



32 7T2 ^ [\ n (fi rj n i)Y 



l + O 



1 



(4.47) 



This leads to the following estimate for the spinodal time for weak coupling A 



f' 2 



ln 



32 7T% 



A(l + 2n ) \a \ 2 fi 



In In 



32 vri 



A(l 



2n ) |a | 2 fl 



O In In In 



(4.48) 



The important point is that the dependence on boundary conditions and the initial distribution is solely logarithmic, 
thus we may expect out predictions to be very robust with respect to changes of the initial conditions. In particular, 
the scale factor at this non-linear time scale is given by 
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where we have used eqs. (|4.3ll and (|4.40|) . 

The amplitude of the long-wavelength modes at the non-linear time, roughly speaking at the end of the phase 
transition is approximately 



\fq(Vnl)\ = 



32 7T2 



A (1 + 2 n Q ) fl 



In 



32 7T2 



_, 1 £1 

4 2i: 



A (l + 2n ) 



(4.50) 



As we will discuss in detail below this non-perturbative scale will ultimately determine the strength of the magnetic 
fields generated during the phase transition. 

During the intermediate time regime the equal times correlation function is approximately 



(*Uv)<S> q , b (v))=6 a . b \a q fe» 2 v 2 e "(# +1 ) ^ 



and its Fourier transform for long wavelenghts is of the form 



(<f> a (x,T])<f> b (0,r))) =S a . b \a \ 2 — /5 3 



In fir] 



(4.51) 



(4.52) 



which determines the time dependent correlation length of the scalar field, 



fa) = - = 2 ■ 



'In (VT^V) 



(4.53) 



This expression is valid in the intermediate time regime rj c < rj < r\ n i during which the non-equilibrium dynamics is 
dominated by the spinodal instabilities. The detailed analysis of the dynamics in refs. |32l l33l l3S| and the discussion 
of the main features presented above can be summarized as follows: 



/i yhi 1/A the mode functions grow exponentially for modes in the 



• At intermediate times fj,^ 1 <C rj < rj n i 
spinodally unstable band q < M(rj). The phase of these mode functions freezes, namely, becomes independent 
of time and slowly varying with momentum. 

• At a time scale determined by the spinodal time the back-reaction shuts off the instabilities and the phase 
transition is almost complete. This can be understood from the following: the backreaction becomes comparable 
with the tree-level term (for 7/ > rjn) when /j 4 ?y 2 . This relation determines that the mean square 
root fluctuation of the scalar field probes the minima of the tree level potential. 

• During the spinodal stage the correlation length of the scalar field grows in time and is given by ea. l|4.53[) . This 
is interpreted as the formation of correlated domains that grow in time, and is the hallmark of the process of 
phase separation and ordering. This correlation length will be important in the analysis of the correlation of 
magnetic fields later. 

• The large fluctuations associated with the growth of spinodally unstable modes of the charged fields will lead to 
current fluctuations which in turn will lead to the generation magnetic fields. Thus the most important aspect 
of the non-equilibrium dynamics of the charged fields during the phase transition is that large fluctuations of 
the charged fields associated with the spinodal instabilities will lead to the generation of magnetic fields. Since 
the modes with longer wavelength are the most unstable the magnetic field generated through the process of 
phase separation will be of long wavelength. Furthermore we expect that the magnetic field generated by these 
non-equilibrium processes will be correlated on length scales of the same order as that of the charged field 
above. 



2. Scaling stage 

A remarkable result of the evolution in the asymptotic regime (when the effective temperature has vanished) found 
in ref. [2(| is that there is a very precise cancellation between the tree level term —fi 2 C 2 (r]) and the back reaction 
XE(i]) in the equations of motion (|4.38|) . The self-consistency condition requires that for a radiation dominated 
cosmology [26) 




(4.54) 
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In this asymptotic regime the solutions of the equations of motion 



are given by 



dr\ 2 



fk(v) = Vv 
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4ry 2 



A k '^^+B k k 2 N 2 (k v ) 



This solution can be written in terms of the scaling variable 

x = k rj 

in a more illuminating form 



f f \ A 5/2 Mx) , „ X 2 N 2 (x) 

fk{v) = A k r) I — g- + B k ^ 3/2 , 



(4.55) 



(4.56) 



(4.57) 



(4.58) 



As discussed in detail in ref. 26], the relevant integrals are dominated by x ~ 1, namely by modes with wavelength 
of the order of the Hubble radius, thus the second contribution proportional to N 2 (x) can be safely neglected at long 
times. 

For x < 1 in the long time regime we can further approximate A k ~ Aq and the asymptotic solution during this 
stage is of the scaling form 



MV)=A 



(4.59) 



Since for x < 3 and large time the modes with small wavevector have the largest amplitudes, these dominate the 
backreaction. The very precise cancellation 14.54fl leads to the following sum rule [2(| 



A 
8^ 



k 2 dk\f k (r,)\ 2 (l + 2n ; 



= A* V 



(4.60) 



Since for large 77 the integral is dominated by soft modes k ~ — — > the distribution function can be approximated 
by n and the amplitude by |^4o| 2 - The sum rule ea. l|4.60[l then leads to the identity 

30?r 3 



\M (l + 2n ) = 



A 



,, 2 TT 2 



(4.61) 



where we used the integral |4£ 



X z 157T 



which is dominated by x < 3. This is a remarkable result: the product |^4o| 2 [1 + 2 no] in the scaling regime does 
not depend on the initial conditions on the evolution, namely it is universal in the sense that it is independent of the 
previous history through the phase transition. This is an important result which will play an important role in the 
power spectrum of the magnetic fields. 

An important consequence of this scaling solution is that the equal time two-point correlation function of the scalar 
field is given by 



($ a (x,ry)$ 6 (0,ry)) = S a>b D(z) ; z 



\x\ 
2 V 



(4.62) 



which reveals that the correlation length is given by the size of the causal horizon [261 . The dynamical evolution during 
the scaling stage is precisely determined by the growth of horizon-sized domains 26] . 

We summarize below the important features of the solutions in the scaling regime that will be used in the compu- 
tation of the power spectrum of the magnetic fields. 
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• For i] >> r] n i a scaling regime emerges in which the mode functions are given by eq. (|4.59() with x — k r\. 
This scaling solution describes the relaxation of long- wavelength fluctuations of the charged fields. Again the 
phase of these modes freezes namely is independent of time. This is important because this fact will entail that 
the retarded self-energy of the transverse photon polarization tensor will give a subleading contribution to the 
generation of magnetic fields. 

• The sum rule ea. l|4.60[) constrains the product of the amplitude times the occupation of the long wavelength 
scaling modes to be given by ea. Q4.6ip . 

• The scaling solution described above is akin to that found in classical models of formation of topological 
defects j^H wl The scaling regime describes the evolution of long- wavelength fluctuations and the adjustment 
of the spatial correlation length of the scalar field to the Hubble radius |26j|. 



C. Gauge field dynamics 

In a high temperature plasma a very important aspect that must be taken into account in the dynamics of gauge 
fields is the electric conductivity, which leads to dissipative processes. As discussed in [27j, the electric conductivity 
severely hinders magnetogenesis, and also introduces the diffusion length scale which could limit the correlation of 
the magnetic fields that are generated. 

In Minkowski space in equilibrium the conductivity is obtained from the imaginary part of the photon polarization 
and it is dominated by charged particles of momenta p ~ T in the loop with exchange of photons of momenta 
eT < k <C T j3|| |37j . A careful analysis including Debye (electric) and dynamical (magnetic) screening via Landau 
damping leads to the conclusion that the conductivity is given by j3(| |37j 

CNT , t . 
a = (4.63) 

with N the number of charged fields and C ~ 0(1). 

In an expanding cosmology an in particular during phase transitions, a more precise assessment of the contributions 
and meaning of the conductivity must be provided. As it was discussed in ref.|23|, the fluctuations of the charged 
fields during a phase transition will have very different behavior if the typical wavevector of these modes is of the 
order of or smaller than the symmetry breaking scale or much larger than this scale. 

Short wavelength modes, those with typical wave vectors much larger than the symmetry breaking scale are insen- 
sitive to the phase transition and are always in local thermodynamic equilibrium (LTE). For short wavelength modes 
deep inside the horizon, the mode functions are of the free field type f q {rf) ~ e rqr) / y/q. 

Long wavelength modes, those with wavectors of the order of or smaller than the symmetry breaking scale undergo 
critical slowing down and fall out of equilibrium during the phase transition. These modes become spino dally unstable 
during the early stages of the transition as summarized above and analyzed in detail in refs. p6l I32L l33l l38(. 

Thus the contributions from the charged particle fluctuations to the photon polarization must be separated into two 
very different regimes: a) the hard momenta p 3> | Af (77) | correspond to charge fluctuations that are always in local 
thermodynamic equilibrium, b) the soft momenta p <C \M(rj)\ fall out of equilibrium and undergo long- wavelength 
spinodal instabilities and enter the scaling regime. 

The contribution from hard momenta will lead to a large equilibrium conductivity in the medium, while the 
contribution to the polarization from soft momenta will contain all the non-equilibrium dynamics that lead to the 
generation of electromagnetic field fluctuations. 

As the instabilities during the phase transition develop, the fluctuations of the charged fields will generate non- 
equilibrium fluctuations in the long-wavelength components of the electric and magnetic fields with the ensuing 
generation of long- wavelength magnetic fields. However the large conductivity of the medium will hinder the generation 
of electromagnetic fluctuations, hence the conductivity must be fully taken into account to assess the spectrum of the 
magnetic and electric fields generated during the non-equilibrium stage[27|. 

We are interested in the generation of long wavelength magnetic fields, namely k « Tr but also k « a 2 Tu, 
since within the astrophysical application, the wavelength of interest for magnetic fields are of galactic scale, while 
Tji corresponds to a wavelength at the peak of the CMB which today is X cm b ~ 1 cm. Thus the physical situation 
corresponds to studying the photon polarization tensor for long-wavelength photons. 

In equilibrium the long-wavelength and low frequency limit (fc, uj — > 0) of the spatial and temporal Fourier transform 
of the transverse polarization is given by 



Hr(fc, uj) = i uj a 



(4.64) 



16 



Thus we write for the full transverse polarization for long-wavelength electromagnetic fields 



II T (?7, rf , k) = o-^6(ij - v') + II non eg(?7, rf, k) 



(4.65) 



The first term above includes the contribution from the hard momentum modes p ~ Tr in the transverse polarization, 
while H noneq (jj,rj' , k) is the contribution from the long wavelength modes which are unstable in the spinodal stage 
and take the scaling form in the scaling regime. Thus in a very well defined sense, the polarization (|4.65|l describes 
the effective low energy theory for the transverse photon field. 

Our strategy is to obtain the non-equilibrium contribution to the spectrum of electromagnetic fields to lowest order 
in a but treating the conductivity exactly. 

In a cosmological space-time, the temperature scales with the inverse of the conformal factor (|4.1U|) and therefore the 
conductivity a = o(jj) becomes time-dependent. If we are interested in time scales where the number of ultrarelativistic 
charge carriers does not change significantly, which is the case that we will consider in what follows, then the time 
evolution of the conductivity is purely kinematic H4.19f) . 

An important effect of the conductivity, as discussed in |2?| , is the introduction of a diffusion scale in the transverse 
photon propagator. The long-time behavior of the propagators for the transverse gauge fields: retarded (R), advanced 
(A), symmetric (H) in absence of non-equilibrium contributions and for ultrasoft momenta k <C o?T 



V^lUr), rf, k) = P«(p) V%>( V , rf, k) 



cbey[T2 ( S ee[23 for details) 
d 2 



^ + k 2 + a( V ) C(r,)j- 
drj z ai] 

dr/ z arj 



V%\n,Tf,k) = 8<j,-rf) ; ©£W) 
V%\r 1 , V ',k) = 5(r,- V ') ; V%\t,,rf) 



= for rj < 7] 
= for rj > rj' 



with the transverse projector 



P'P J 



(4.66) 



(4.67) 



Due to eq. Ij4.19|) the comoving conductivity or = <r(rj) C{rj) is an invariant quantity in the regime in which the 
number of ultrarelativistic charge carriers is constant. The estimate given by ea. (|4.22|) clearly indicates that during 
the radiation dominated era between reheating and the QCD phase transition, or rj 3> 1. 

Then for k <C or (which is certainly fulfilled since the relevant wavevectors are i;«r« T/a ~ or) and rj 3> 1/or 
we can safely neglect the second order time derivatives in eos. (l4.66|l . leading to the following equations, 

(4.68) 

VP(r,,rj';k) = ^- . (4.69) 



2>£> (r,, r,',k) = V$ (77, rj'; k) % - rf) , (77, rf, k) = V { ° ] (rj, rf; k) 9{rj> n) 



V%\ri,rf,k)=i 



or 



or 



V. MAGNETIC FIELD SPECTRUM 

As discussed in detail in reference [22}, the quantity of astrophysical relevance is the correlation function 

< ^(77,^)^(77,0) > p , (5.1) 

where the sum on repeated indices is understood. B(rj,x) above is a Heisenberg operator and the expectation value 
is in the initial density matrix. From this quantity, the spectrum of the magnetic field is obtained in the coincidence 
limit 

S B (v,k) = \ km / d 3 x < {B\rj,x),B\rj'M > p e^ s , (5.2) 



17 



where { , } denotes the anti commutator. And from $#(77, A;) we can extract the physical magnetic energy density 
stored on comoving length scales larger than a given L 



A PB (L, V) = ^2 I " k 2 S B (m k) dk . (5.3) 



where we have restored the powers of the scale factor arising from the transformation to conformal time. Denoting by 
Apb(L,t]) the contribution from the non-equilibrium generation (subtracting the local thermodynamic equilibrium 
contribution) , a quantity of cosmological relevance to assess the relative strength of the generated magnetic field is 
given by the ratio of the power on scales larger than L to the energy density in the radiation background 

r(L,v) = ^^, (5.4) 
Pivl) 

where 

7T 2 T 4 

P, = ~jf (5-5) 

is the comoving energy density in the thermal equilibrium background of photons. 

The physical energy densities ApB, P hys(L, rf) , p ltP h, ys are obtained from the comoving expressions above by rescaling 
p — » p/C 4 (rj) as can be seen from the conformal rescaling <|3.8fl . Thus the ratio r(L,rj) would be a constant in the 
absence of non-equilibrium generation or dissipative processes. Hence the time dependence of the ratio (15.411 only is 
solely a consequence of the non-equilibrium generation mechanisms or dissipative processes (such as magnetic diffusion 
in a conducting plasma) but not through the cosmological expansion. 



A. Exact formulation 

In terms of the transverse component of the gauge field A l T (rj, k) — V %3 (k) A 3 (r/, k) with the transverse projection 
tensor = 5y — ^l 1 the spectrum S B (r], k) can be written as (with an implicit sum over indices) 

Ss(v,k) J d 3 xe lts k 2 V%{r,^-x)\ n=ii , (5.6) 

in terms of the symmetric correlator of the transverse gauge field 

V%( v ,r,i;x) = V i3 >(r 1 ,r 1 ';x) +V i3 <(r ? ,r ? ';x) 
^>(ri,Tf;g) = i<A i T {r,,x)Aj r {r 1 ',$)> ; X>« <(r,,ri';x) = i < 4(t/,0)^M > (5-7) 

Where the correlation functions in (|5.7fl are Wightmann functions without time ordering computed in the initial 
density matrix. The main reason for introducing the definitions above is that there is a well established framework 
for obtaining these correlation functions in non-equilibrium quantum field theory as discussed below. 
^From Sb{i1i k) we can extract the magnetic energy density 

/d 3 k 
-—-r^S B (ri,k), (5.8) 

where the ultraviolet behavior is understood to be regulated in some gauge invariant manner, for example with 
dimensional regularization. 

^From the magnetic energy density we can define an effective magnetic field B e ff(r]) such that 

^B 2 eff ( V )= P B(v) (5.9) 

In the non-interacting, thermal equilibrium case, which will be relevant to compare the energy density in the generated 
magnetic field to that in the thermal radiation background, 

S { °\k) =k [l + 2n(fc)] (5.10) 
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where n(k) is the Bose-Einstein distribution function. Subtracting the irrelevant vacuum contribution, the effective 
magnetic field is then given by the Stefan-Boltzmann form 

*eV 0) = ^ ■ (5-11) 

The phase transition generates a dynamic effective magnetic field A_B e //(?y) through the interaction between the 
charged fields and the electromagnetic field. Hence a priori we would expect that AB e //(r;) can be obtained system- 
atically in a power series expansion in the electromagnetic coupling a, namely AB e ff(i]) = aAB < e 1 jj : (rj)+a 2 AB < fjj.-(rj) + 

This would be the case were it not for the fact that the DC conductivity is important to estimate reliably the 
amplitude and the correlations of the generated magnetic field. As it will be discussed below in a high temperature 
plasma the conductivity is determined by fluctuations of the charged fields with typical momenta p ~ T. These short 
wavelength modes are in LTE while only low momentum modes of the charged fields will fall out of equilibrium and 
undergo spinodal instabilities during the phase transition. The DC conductivity is nonperturbative and non-analytic 
in a and to leading logarithmic order it is approximately given bv |36ll37| 

CNT (5.12) 



a In -Kt 

aN 



with TV the number of charged fields and C a numerical constant of order one|36l l37|. 

Thus while we cannot provide an expression for the generated magnetic field as a power series expansion in a 
because the presence of the conductivity prevents such expansion, the strategy that we pursue here is to treat 
the long-wavelength non-equilibrium current fluctuations in perturbation theory in a while the short wavelength 
contribution will be accounted for in the conductivity. This will be discussed below in detail, but the important point 
of this discussion is that the long-wavelength fluctuations that lead to the generation of the magnetic field will be 
treated to lowest order in a and to leading order in the large N limit. 

The reliability of this expansion will be guaranteed if 

<«»> 

In cosmology the important quantity is the energy density in long-wavelength magnetic fields on scales L equal to or 
larger than the galactic scales. Assuming rotational invariance, we introduce the energy density of the magnetic field 
generated by the non-equilibrium fluctuations 

A PB (L) = — k 2 AS B (k,t) dk . (5.14) 

2tH Jo 

The quantity of cosmological relevance is 

Pi 

where p~ ( = 7r 2 T 4 /15 is the energy density in the thermal equilibrium background of photons. 

As discussed in section IV, the physical magnetic field in a cosmological background is diluted by the expansion as 
B 1/C 2 (n) with C{rj) the scale factor. Therefore, in the absence of processes that generate or dissipate magnetic 
fields the ratio (|5.15|l would be constant under cosmological expansion because of the redshift of the temperature T ~ 
1/(7(77). I n a cosmological background the ratio (|5.15|l will only depend on time through the generation or dissipative 
mechanisms (such as magnetic diffusion in a conducting plasma) but not through the cosmological expansion. 

For r > 10 -34 the linear (kinetic) dynamo may be sufficient to amplify cosmological seed magnetic fields and for 
r > 10~ 8 the co llap se of protogalaxies with constant magnetic flux may be sufficient to amplify the seed magnetic 

fieidsnssmiii. 

In the cosmological setting this ratio is approximately constant if the conductivity in the plasma is very large after 
the non-equilibrium dynamics has taken place. The approximations invoked to estimate the spectrum of the generated 
magnetic field will be reliable provided r(L) <C 1. 
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1. N on- equilibrium dynamics of electromagnetic fluctuations 

The formulation of quantum field theory out of equilibrium is by now well established in the literature and we 
refer the reader to refs [3^, H3 ED for details. The generating functional of real-time non-equilibrium correlation 
functions can be written as a path integral along a contour in (complex) time. The forward and backward branches 
of this contour represent the time evolution forward and backward as befits the time evolution of an initial density 
matrixj32,[3!|. Gathering all bosonicfscalar and vector) fields generically in a multiplct the generating functional 
of correlation functions is given by [89L EflL El) 

Z(J + , J_) =y"[P*+P*_] e */«J**W*+)-A*-)+J+*+-J'-*-} . (5.16) 

where functional derivatives with respect to the sources J± lead to the non-equilibrium real time correlation functions. 
The doubling of fields with labels ± is a consequence of the fact that the non-equilibrium generating functional 
corresponds to forward (+) and backward (— ) time evolution and suggests introducing a compact notation in terms 
of the doublet (\&+, and the metric in internal space 

<** = (J \) = tat » (5-17) 

where the labels a, b = ±. This notation is particularly useful to obtain the non-equilibrium version of the Schwinger- 
Dyson equations for the propagators |39j . 

In the case under consideration the main ingredients for the non-equilibrium description are the following 

• Transverse photon propagators 

The real time Green's functions for transverse photons are given by 

(4^(f,t)^j(P,0) = -ij <(M>0 e-^-*') , (5.18) 

where the explicit form of D?- (fc; t,t') is 

V±+(k;t,t')=V ij (k) [P>(A;;t,t')e(7 ? -t')+I?<(fc;t,i , )0(^-i)] , (5-19) 
V7.-(k;t,t')=V ij (k) [V^t^eir}' -t) +Vf j (k;t,t')e(r, -t')] , (5.20) 
V±-(k;t,t') = Vij{k) V<(k;t,t') ; Vr+(k-t 1 t l ) = V^k) V>(k;t,t') (5.21) 

and Vij(k) is the transverse projection operator l|4.67[) 

• Scalar propagators: 

(* r W(S,t)*W@,ir)) = ~i5 rs J ^ Gt{r,,t') e-^< s -^ , (5.22) 



G+ + ( V ,t') = G>( V ,t')e( V - t') + G<(r),t')G(r]> - t) , (5.23) 

G--( V ,t') = G>(vXMrf-t) + G<(r ] ,ti)e(r ] ~t'), (5.24) 

G+-( Vl t') = G<(v,t') ; G^ + (r,,t') = G>(r,,t>), (5.25) 

G£(t?,0 = |{[l + n fc ]/ fe (7 ? )/ fe *( ?7 / )+n fe / fe (V)/ fc *W} , (5.26) 

0^(77, = | {[1 + n k ] M V ')r k ( V ) + n k f k (rj)r k (v')} ■ (5.27) 



For the scalar propagators we have used the expansion in terms of mode functions given by eq. (|4.29|1 and the 
expectation values given by eq. 14.32JI . 

The scalar propagators given above imply a non-perturbative sum of cactus-type diagrams when the expectation value 
of the a field vanishes. These propagators are depicted in fig. 
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FIG. 1: Scalar propagators in leading order in the large N. For vanishing expectation value of a these are given by a sum of 
cactus-type diagrams. 



2. N on- equilibrium Schwinger-Dyson equations 

In this section we derive an exact expression for the spectrum of the magnetic field from the non-equilibrium set 
of Schwinger-Dyson equations. 

The full non-equilibrium propagator for the photon field is obtained from the non-equilibrium effective action 
resulting from integrating out the charged fields, and which up to quadratic order in the photon field is given by 

r[A±] = J d'xd'y {Ay(x) [d^S\x - y)c ab S l3 + IL% T (x, y)] A b . T (y)} (5.28) 

with an implicit sum over all repeated indices. 

To lowest order in a and to leading order in the large N expansion, the non-equilibrium contribution from the 
scalar fields to the photon polarization is given by the two diagrams shown in fig. |2 




(b) 



FIG. 2: Photon polarization to lowest order in a and to leading order in the large N. The scalar loops are in terms of the 
scalar propagators displayed in fig0 

The tadpole diagram (a) in figEl gives a contribution that is local in time, thus we write the polarization 

n$ tT (x,y) = n tad ( v ) S i (x-y)+U^ T (x-y,x°,y ) , (5.29) 

where the form of the local (tadpole) contribution Ii tad (rf) in terms of the metric c ab is necessary for consistency 
with the form of the effective action since the kinetic term is proportional to the metric 39] . The contribution II is 
non-local, and in equilibrium it features absorptive parts. The different components a, b — ± arc defined in the same 
manner as the Green's functions i|5.19|l - i|5.21[l an d l|5.23|) - (|5.25|) . The contribution to order a to the non-local part of 
the polarization is depicted in diagram (b) in figEl 
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The photon propagator given in eq. I|5.18jl is the inverse of the operator in the quadratic effective action above and 
obeys the non-equilibrium version of the Schwinger-Dyson equation |39| 



d 4 y [d^5\x - y) c ab 5 l3 + U$(x, y)] V{ k c { Vl z) = 5 k c a c S\x - z) 



(5.30) 



From the expressions H5.19M5.21j l f° r the different components of the photon propagator it proves convenient to 
introduce the following combinations 



Vf j (x,y)=Vf j (x,y)-Vf j (y,x) 
V? j (x ) y)=VC j (x,y)e(x°-y ) 
Vf j (x ) y)=V? j (x,y)e(y -x ) 
T£(x,y)=V>(x,y)+V<(y,x) 



antisymmetric propagator 
retarded propagator 
advanced propagator 

symmetric propagator. 



(5.31) 
(5.32) 
(5.33) 
(5.34) 



and similar definitions for the polarization II. 

Using spatial translational invariance we can simplify the above form of the Schwinger-Dyson equations by intro- 
ducing the spatial Fourier transform of the transverse propagators 



VtU V ,7 1 ';x-x') 



d 3 k 
(2^)3 



and introducing the spatial Fourier transform of the transverse polarization tensors defined by 



d 3 k 



(5.35) 



(5.36) 



We can now write down the Schwinger-Dyson equations for the transverse photon propagators. Since there are only 
two independent functions D > < the Schwinger-Dyson equations Ij5.3()jl can be simplified to a set of two equations 
for a particular combinations of these. Some straightforward algebra leads to the following set of non-equilibrium 
Schwinger-Dyson equations (we have suppressed the subscript T but the equations below refer to the transverse 
components) 



— + e + Il tad ( v ) ) 5(ri-Tn) + iL R (v,Th;k) 



and 

with the definition 



dr)i 



(5.37) 



fff^i; 7 ?'; fc ) = / d-qi Tl H (77,771 ;k)V A (rii, 77'; k) . (5.38) 



U R (r),rn;k) = ^(rj^k) -H<{ri,rn;k) 6(771-*) 



(5.39) 



A remarkable aspect of this set of equations is that the right hand side of eq. Ij5.38jl . namely the inhomogeneity in 
the set of equations, only involves the non-local contribution to the polarization, which to lowest order in a is given 
by diagram (b) in fig. [21 This is a result of the form of the local (tadpole) contribution which is proportional to the 
metric c ab . This is an important point to which we will come back later in the discussion. 

The antisymmetric propagator T>c{n, 77"; k) is the odd solution of the homogeneous evolution equation 



+ ^ + ntad ^) Vc ^ V"; k) + ^ dt' U R ( V , 77', k) V c (rf, 77", k) = , 



with the constraint 



drj 



= 1 



(5.40) 



(5.41) 



This relation guarantees the correct equal time canonical commutators. 
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Since the kernel of the integral equation on the left hand side of eq. <|5.38[1 is the same as for the equation that 
defines the retarded propagator l|5.37[) . the solution to eq. 1)5. 38[) is given by 

V H (ri,rj';k) = J dm Jdrj 2 V R (rj,^i;k)n H (r]i,i]2;k)V A (r]2,r]';k) + F(r],ri';k) (5-42) 

where D B (m t\'\ k) is the solution of eq. (|5.37() and the function F(m rj'; k) symmetric in r\ and r( is a general solution 
of the homogeneous equation 



dm 



dt\ 



+ k 2 +n tad (ry) 5(r 7 -7 7l )+n fl (r ? ,r 7l ;fc) 



F( Vl ,r]',k) = 



(5.43) 



The homogeneous solution F can be constructed systematically and its physical significance will be discussed below. 

Now are now in position to provide the final expression for the spectrum of the magnetic field. From the expression 
(15.6)1 and the solution found above, we find 



S B (r],k) = -ik 2 V H { v ,ri';k)\ 



7J = 7j' 



(5.44) 



The expression (|5.42|l can be simplified further by taking into account the theta functions in the definitions of the 
retarded and the advanced (T>r and T) a) propagators [see eqs. 1)5. 32)1 - 1)5.33)) ] as well as the antisymmetry T>c(m rj'; k) — 
—T>c{i]', t; k) of T>c leading to 

S B (v,k) = i k 2 I ( drji f dr} 2 Vc(v,V2;k)Vc(mVilk)TL H (m>W,k)-F(mt;k)\. (5.45) 

Since the product T>c{mm\kYDc(r\,m\k) is symmetric in the exchange m «-> m we can replace TLH(r]i,T)2;k) by 
2n>(?7i, 772; k) and write the final form of the spectrum separating, for further convenience, the contribution from the 
inhomogcneous and homogeneous solutions to Sb- 



S B ( V ,k) = S B (r ) ,k) + S B I (r 1 ,k) 

S B (r),k) = 2i k 2 / drji / dr] 2 £>c(»7, V2; k) V c (r], 771; k) n>(r/i, 772; k) 

Jr>R JriR 



S§(r,,k)=-ik 2 F(rj,t;k) 



(5.46) 



where r/ji is some initial time before the phase transition and n> does not include the local tadpole contributions, it 
is the non-local part of the polarization. 

There is an important aspect associated with the homogeneous solution F(r),r)',k) and its contribution to the 
spectrum of generated magnetic fields S^(n, k). This aspect is revealed by noticing that the expectation value of the 
transverse gauge field -4t(?7, k) = (At{ti, fc)), namely the mean field obeys the same homogeneous equation of motion 
as F{r},r}',k), 



dm 



dt 2 



+ k 2 + Il tad (r,)) 5(r,-m)+^R(v,m;k) 



A T (v,k) = 



(5.47) 



(we have suppressed the vector indices to avoid cluttering of notation) . 

Thus the homogeneous solution F(r],r)',k) can be constructed out of the independent solutions of the mean field 
equations of motion l|5.47|l . The main reason that we bring up this point is to highlight that the solutions to the 
mean-field equations of motion are only part of the contributions to the generation of magnetic fields through non- 
equilibrium processes. However, as it will be discussed in detail in the next section, this contribution can be neglected in 
the present case of non-equilibrium spinodal decomposition in many circumstances, and the term S B (n, k) dominates 
for late times. 



3. Electric fields: 

For completeness we now address the generation of electric fields. The importance of generation of electric fields is 
mainly related to the question of equipartition. It is often assumed that the energy density stored in electromagnetic 
fields is equally partitioned between electric and magnetic field components, namely between temporal and spatial 
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gradients. While this is usually the situation in equilibrium, it is not necessarily the case strongly out of equilibrium 
such as the situations envisaged in this article. 

The electric field is the Hamiltonian conjugate field to the vector potential, and its transverse component is given 
by E l T (r\, k) = —A l T (rj, k), therefore the equal time correlation function of the electric field is given by 



S Et (v,x) =< Ai r {7 1 ,x)A l T {i 1 ,6> c = d v d, r < i{A^(»7,£),A^(j/,0)} > t 



(5.48) 



where just as in the case of the magnetic field we wrote the equal time correlator as the symmetrized connected 
two-point correlation function. 

Following the steps leading to eq. 1(5.45(1 for the magnetic field, and using the following identities 

drt'^BiViVi'ik) = d v T> c (v,m; k ) Q{n - Vi), c^XUfe?/; k) = -d v iT>c{rf k) 0(i] 2 - rf) 

we now find 

S ET (-n,k) = d n d n , hi ['dm [ DciVtmifyDciv^Viltyn^^^k) - F(r),ri',k)\ (5.49) 

The number of photons produced through the non-equilibrium processes is given by 

= S ET (ri,k) + S B ( V ,k) 1 
2k 

after summing over the two polarization states. 

The final form for the spectrum of magnetic and electric fields given by eq. 1)5.46(1 - 115.49(1 as well as the number of 
photons given by eq. 15.50(1 are the main tool to compute the spectrum of electromagnetic fields generated during 
non-equilibrium processes and one of the main results of this article. 

We emphasize that these expressions are exact and general and apply in the cosmological setting. They allow to 
study the problem of the generation of magnetic fields through non-equilibrium processes in general situations. 

Furthermore, the final expressions for the spectrum of electromagnetic fluctuations ((5.461 l5~4*9"|l is valid more gener- 
ally in spinor electrodynamics since it only involves the full polarizations and the Schwinger-Dyson equations for the 
correlation functions. 



4- Spectrum of fluctuations in equilibrium 



Before focusing on the study of the magnetic field generated during non-equilibrium phase transitions, it is both 
illuminating as well as important as a consistency check to address the case of thermal equilibrium in Minkowski 
space-time. In this case the simplest manner to compute the spectrum is by usin g th e imaginary time or Matsubara 
formulation, where the Euclidean transverse photon propagator is written as |42l l43j 

/d 3 k - -» 
« e ik-g T ^p«(fc)p T ( Wn ,fc) e -^n(T-0 . Wn = 2 tfn (5.51) 



with 



V T (u n ,k)= fdko^^. (5.52) 

J k ~ lUJ n 

and pT{k,ko) is the spectral density for transverse photons which is an odd function of kg. The spectrum for the 
magnetic field is obtained from the equal time limit of the Euclidean propagator, and is therefore given by 

S B (fc) = 2fc 2 T^2Mw„,fc) . (5.53) 

nez 

The sum over the Matsubara frequencies can be performed using the methods described in p2L |43| leading to 



S B (k) = 2 dujk z n(w) p T (u, k) . (5.54) 
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The general form of the spectral density in terms of the transverse polarization is given by (42, 4. 

fl , \ 1 Imn T (fc,fc ) 

PT{k,ko) = 5 =■ (5.55) 

TT [fc 2 - fc 2 _ R e n T (fc, fc )f + [ImIT T (fc, k )] 2 

The spatial and temporal Fourier transform of the retarded and advanced transverse photon propagators are given 

by mm 



' '■' • • ' ' - ' -^ 2 + fc 2 + RenT(fc,cj) + iImn T (fc,w) 

1 

-cj 2 + fc 2 + ReIT T (fc,cj) - i ImII r (fc,cj) 



W )=P tf (fc)2> A (fc,w) ; V A (k,u)= 2 ■ ;2 ; ^_ (u - , TmTT /; _ A (5.56) 



We are now in conditions to establish contact with the non-equilibrium result of the first section. In equilibrium the 
polarization and the propagators are functions of the time differences and we can then take the Fourier transform in 
time of T>h{t], r)' , k) given by eq. 15.42fl 

T> H (k,co) =V R (k,Lo) p>(i:,w)+n < (J:,w)]P4fc,w)+f(fc,w) (5.57) 

where we have explicitly written TIh in terms of II<,>. In equilibrium the detailed balance (or KMS) condition 
relates these components of the polarization to the imaginary part |42l |4^| 



n>(fc,w) = -\l+n(u})] ImII(k,u) ; n<(&, w) = - n(ui) lmU{k, u) (5.58) 

with n(u>) being the Bose-Einstein distribution function. 

Furthermore, the Fourier transform of the homogeneous solution F(r),rj',k) obeys the Fourier transform of eq. 
(|5.43|l namely 

V^iuj.k) F{uj,k) = (5.59) 



Combining all the above ingredients and using the fact that ImII(A;,a;) is an odd function of o;|42j,|43| we are led to 
the conclusion that eq. H5.44|l for the spectrum becomes 

S B (k) = [ dw {2 k 2 rc(cj) p T (u>, k) + k 2 F(u), k)} (5.60) 



with Pt(w, k) given by eq. I|5.55|l . 

This result differs from that obtained via the equilibrium propagator eq. Q5.54J1 by the contribution of the homoge- 
neous solution F. 

In order to understand the source of the discrepancy between the two formulations, namely the contribution of the 
homogeneous solution F(k,u), let us focus on the defining equation for F(k,u>) l|5.59|l . 

In order for a non-vanishing solution to this equation to exist, from the expression for Vn(k,to) we infer that 

uj 2 - k 2 -ReU T (k,uj) = (5.61) 
ImII T (fe,w) = (5.62) 

Eq. (|5.61|) determines the dispersion relation of quasiparticles and eq. Ij5.62|l determines that these quasiparticles 
must have zero width, i.e, the solutions of the homogeneous equations are propagating quasiparticles with zero width 
and a dispersion relation given by (|5.61(l . The homogeneous solution is therefore 

F(k,u) (x 5{lo 2 - k 2 -ReU T (k,uj)) (5.63) 

In a non-perturbative resummation of the Dyson series for the propagators and the spectral density, the limit when 
ImllV — ► 0^ leads to 

p T -> tS(lo 2 -k 2 - Ren T (fc, uj)) (5.64) 

and we recognize in this case that the possible contributions of the homogeneous solutions are accounted for in the 
expressions l|5.55|l and arise from the propagating states, and to zeroth order account for the contribution to the 
spectrum of the magnetic field from free photons. 
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However, in a perturbative expansion perturbation theory is only reliable away from the single particle poles in the 
propagator and the limit Imlly — ► 0^ will miss the contribution from the isolated poles. The homogeneous solution 
F(k, u) gives the required contribution, thus guaranteeing the consistency of the perturbative expansion. Thus while 
the homogeneous solution F must be included in a perturbative calculation of the spectrum, it will be accounted for 
in a non-perturbative computation that includes a Schwinger-Dyson resummation for the full propagator and should 
not be included in the computation of the spectrum. 

In a plasma the true degrees of freedom are quasiparticles, in particular for the power spectrum of magnetic 
and transverse electric fields, the important degrees of freedom are transverse plasmons. In the hard thermal loop 
approximation)^ |43| the transverse plasmons contribute to the spectral density pr(k,uj) a pole term of the form 
Zx(k)[6(u} — u> p (k)) — S(lj + u> p (k)] with u> p (k) being the plasmon dispersion relation |42l l43| with lj p (0) = eT/3. For 
k ^> eT, ui p (k) ~ k ; Zx ~ I /2k and we recover the power spectrum of free field theory. However for k < eT the power 
spectrum of the magnetic field reveals the presence of collective degrees of freedom, Ss(fc) ~ kZr(k)[l + 2 n(u> p (k)] 
for e 2 T <k< eT. 

For ultrasoft frequency and momenta k,uj <C c?T which is the case of relevance for the cosmological case, the 
effective low energy form of the gauge field (transverse) propagator (well below the plasmon pole) is 

pk g) 2 i ■ ( 5 - 65 ) 

U) z — k z — IO10 

leading to the small frequency, long-wavelength form of the spectral density 

1 a to 

7T (lo 2 - k 2 ) 2 + (crw) 

Therefore, for k <C cr and lo <C k 



1 a lo 

p T (k,Lo)~-— — . (5.66) 



k 2 n(Lu) PT (k,Lo)^-^—t ; T= — . (5.67) 

Thus for long wavelengths k <C o ~ T / a the spectral density features a pole at zero frequency leading to a long- 
wavelength power spectrum 

5 B (fc< ( j)=r, (5.68) 

coming from the first term in ea. (|5.60|) . 

The plasmon pole contributes to ^(fc) through the homogeneous solution F in ea. H5.60[) yielding, 

k 2 2k 2 T k 2 

SBplasmon - 7^[1 + 2n(uJ p (0)] — - — (5.69) 

u p (0) Lu z (0) e 2 T 

where we have used the long- wavelength limit of the dispersion relation u p (k) and residue Zx(k)^^, |4^. The 
contribution (|5.69|) is clearly much smaller than l|5.68|l in the long-wavelength limit k -C eT. Hence, the magnetic 
field power spectrum in equilibrium is Ssik) — T. 

Having clarified the role of the homogeneous solution and the form of the power spectrum within the more familiar 
equilibrium setting, we return to the non-equilibrium situation. 



5. Spectrum out of equilibrium 

Using the results obtained above, in particular the final expression for the magnetic field spectrum, eqn. 15.461 we 
are now ready to study the spectrum of the magnetic field generated by the non-equilibrium evolution of the charged 
scalar field. 

While the final result for the spectrum given by eqn. (|5.46|) is exact, in order to make progress and obtain an 
estimate for magnetic and electric field generation we have to make approximations. In what follows we will obtain 
the spectra for magnetic fields to lowest order in a and to leading order in large N. 

To this order the tadpole (local) and bubble (non-local) contribution to the polarization are given respectively by 
diagrams (a) and (b) in fig. [21 Their explicit expressions are given by 

U tad (t) = ^e 2 N J ^ G> (t,t,q), 
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nij,>(j7i, 772,fc) = ie 2 N 
which leads to the transverse components 
n>(f7i, 772,fc) = 2ie 2 N 



d 3 g 
(2nf 



(2<r + ki) (2% + kj) G> (771, 772, g) G>(r ?l! 772, 1? + 



(5.70) 



2 v / rf3 g . 2 , j :.; 

(2tt) 3 



9 2 (1- cos 2 6<) G>(riuri2,q) G>(f?i, 772, W + k\) , 



n<(r?i,?72,fc) = n>(?72, 77i,fc) 



(5.71) 



where cos 9 = q ■ k. The scalar propagator G>(ry,?/) in terms of the mode functions f q (t) that satisfy the evolution 
equations (|4.3U|) is given by [see eq. I|5.26|l ]. 



g> (771,772;*;) = - [(1 + n k )fk(rn)fk(m) + nkfk( T n)fk(v2)] 



(5.72) 



and we assumed that before the phase transition the scalar fields have occupation n^. Inserting this expression in 
Sg(rj, k) in eqn. i|5.46[l we find 



S I B {r ] ,k) = e 2 N 



d 3 q 
(2tt) s 



g 2 (l-cos 2 6) 



(H-n g )(l+n lr+jf| ) 



+(i + 0»V+£| 



dru fcX> c (77,77i,fc)/ g (77i)/*_. + - | (77i) 



1R 



+ n i n \q+k\ 



d77x k T> c (ri,Tlu k )P q (Vt) f* ?+ ^(Vt) 



(5.73) 



The different terms in this expression have a simple kinetic interpretation. The first term, proportional to (1 + 
n q )(l + Hig._|_a) corresponds to the process of pair creation along with the creation of a photon, the second and third 
terms, proportional to the combinations (1 + n)n correspond to the process of bremsstrahlung, and the last term 
proportional to n n corresponds to pair annihilation with the emission of a photon. These processes are depicted in 
figure El 





(a) 



(b) 



(c) 



FIG. 3: Processes described by S B . (a): pair annihilation, (b): pair production, (c): bremsstrahlung 

These processes are forbidden in equilibrium by energy momentum conservation, however in an expanding cosmology 
and strongly out of equilibrium these processes are not only allowed but they give the leading order contribution to 
magnetogenesis . 

The contribution from the homogeneous solution S B (77, k) is given by 



S^(r 1 ,k)=-ik 2 F(r 1 , m k) , 
where F{rj,r]';k) satisfies the homogeneous differential equation 



drf- 



dr] 



F(r), V ',k)+ d Vl [TL tad (r n )5(r 1 -j n )+Tl R ( V ,r ll )]F(r ]1 ,r 1 ',k) = 



(5.74) 



(5.75) 



with U tad (r]), IIr(77, r)') being the one loop tadpole (local) and retarded (non-local) contributions transverse 
polarization[27|. 

We note that the function F(j], 77'; k) obeys the same equation as the transverse gauge mean field |27||. but as it will 
be argued in detail below, its contribution to the spectrum of primordial magnetic fields generated during the phase 
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transitions is subleading in the scalar coupling constant A. The equation of motion l|5.75[l can be solved systematically 
in an expansion in powers of the non-equilibrium polarization, 

F(n, m k) = (n, m *) + (v, m k) + o(a 2 ) , 

F^(r},tf,k)=V B (r},tf,k) , 

F (1) (l.l'.*)= / d m V c (r,-m,k) [ d m \n tad (r} 2 ) S{ m - m ) + TL R ( m , m - k)] V H (r? 2 - 77', k) + 



nn 



(5.76) 



where Tl tad , IL R are the tadpole (local) and the retarded contribution from the one-loop transverse photon polar- 
ization respectively (for details see[27|1. 
These are given by [2^ 

II tad (?7) = -ie 2 N 

with 

II>(77i,7&,fc) = 2ie 2 N 



d 3 g 
(2tt) 3 



G> fa, ??, g) , Il R (r]i,r] 2 ,k) = rPfai,?^, k) - IT^fai, r? 2 , k) 6(771 - 772) (5.77) 



q 2 (1 - cos 2 6>) G> (771,772, g) G> (771, 772, |<f+fc|) , n<(77i,?72,/c) = n> (772 , 771 , fc) 



The scalar propagator (?>(?;, 77'; fe) is constructed from the mode functions ^(77) that satisfy the mode equations (|4.42|) 
and is given by l|5.72|l : therefore 



d 3 q 
(2ny 



q 2 (1 - cos 2 



Im 



G> (771,772,9) G> (771,772, |g + fc|) 8(771-772) . 



(5.78) 



This expression for the retarded self-energy must be contrasted with that of the contribution from Sg(r),k) which 



requires the real part Re 



G>{m,V2,q) G>(t7i,?72, \q + k\) 



This is an important difference, the long wavelength 



modes of largest amplitude in either phase given by (|4.45|l or by (|4.59|l are such that their phases are frozen, namely 
they do not depend on time, therefore the products f q {'r]i)f\p+^\ (772) with only the growing mode solutions are real 
and such products will contribute only to S^(i], k). This freezing of phases is a consequence of the classicalization of 
the scalar field fluctuations [26| . 

We now argue that the contribution from Sg is subleading. First of all, the term F^ ' (77,77; k) in eq. lj5~7T)|l is the 
solution of the homogeneous equation in absence of non-equilibrium fluctuations and leads to the local thermodynamic 
equilibrium contribution to the power spectrum, which is independent of the non-equilibrium generation mechanisms. 
This contribution has been analyzed in detail in ref.j^ and will be subtracted. In what follows we focus solely on 
the contribution from the non-equilibrium fluctuations. 

For intermediate times after the phase transition during the spinodal stage r\ c < 77 < 77,^ , the long- wavelength mode 
functions are approximately given by ea. 14.45fl . 

Near the end of the phase transition for 77 ~ 77' ~ 77,^ the leading order time dependence of the scalar Green's 
functions is approximately given by 



fk(v)fk{v') °c 



A 



(5.79) 



where we used ea. (|4.50() . Thus, the contribution from the tadpole (local term in the self-energy) is of the order 



n (r))v~vm 00 ~r + subleading . 
A 



(5.80) 



This estimate is consistent with the fact that the tadpole contribution is e 2 < <!><!> > and near the end of the phase 
transition the mean square root fluctuations of the scalar field probe the vacuum state, namely < <f> 2 >^ /J 2 /A. Since 
the phases of these modes are frozen, there is no contribution from the leading order to the retarded polarization, since 
it requires the imaginary part of the product of propagators as displayed in eq. (|5.78(l . Because of this cancellation 
of the leading term, the contribution from the retarded polarization bubble is of the same order as that of the 
tadpole p7Ll46l|. 



U R oc 



(5.81) 
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A similar argument based on the sum rule (|4.60() leads to the same conclusion in the scaling regime. 

The contribution from S B is in both cases of 0(1/ A 2 ) since each long- wavelength mode function is of order 1/y/X 
at the end of the spinodal stage or, by the sum rule l|4.60|l in the scaling regime. Thus we can safely neglect the 
contribution from Sg to the magnetic spectrum. 

The general expression (|5.73[> can be further simplified by noticing that during most of the phase transition the 
mode functions are frozen: therefore the four terms in (|5.73(l can be collected together and the leading contribution 
to the power spectrum generated by non-equilibrium fluctuations can be written as 



S B {r),k) = (l + 2n ) : 



aN k 2 



q i dq d(cos9) (1 



:>ii 



fq(Vl) /| 9 -+fe|(^l) dr K 



(5.82) 



where 9 is the angle between the vectors q and k and where we have replaced 

(l + 2n,)(l + 2n |f+g| )~ (1 + 2n ) 2 , 



since as highlighted in section IjlVB ljl the dynamics during both the spinodal stage as well as the scaling stage is 
dominated by the long-wavelength modes that acquire non-perturbatively large amplitudes. 

The final form of the power spectrum generated by the non-equilibrium dynamics given by eq. I|5.82|l is the basis for 
the study of primordial magnetogenesis during the different stages after the phase transition. 



B. Magnetogenesis during the spinodal stage 

The long- wavelength mode functions in the spinodally unstable band are given by the expression l|4.45|l . 

The integral over r/i for large rj can be computed integrating by parts in ea. H5.82[> as an expansion in I /(fir]) 2 . 
The integral is dominated by the upper limit, which leads to the cancellation of the exponentials that contain the 
conductivity. 

The integrals over momenta and angles in ea. (|5.82|) can be done straightforwardly when the mode functions are 
given by ea. H4.45(l . Thus from ea. l|5.82|l whe obtain the following expression for the spectrum of magnetic fields 
generated by the non-equilibrium fluctuations 



S B (k,r)~ T] n i) 



51271-5 N a k 2 

A 2 o\ A 4 e(vm) 



1 + 



where is given by ea. (|4.53|) . In obtaining this result we used the following 



[1 + 2 n ] 2 \a Q \ 4 



[see ea. (|4.48l) ] and the identities^ 
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(5.83) 



(5.84) 



(5.85) 



and 



/ qdq e~ e q2 [Z 2 q k cosh(f 2 q k) - sinh(^ 2 qk)] = ^- k 3 £ 
Jo ° 



(5. 



Notice that the magnetic field spectrum (|5.83(l is independent on the amplitude |ao| and on the initial occupation 
(1 + 2?i,o) 2 . Therefore this result is quite robust. 

This result is the same as for the Minkowski space-time (see eq.(7.47) in ref.j^), except for a multiplicative factor 
fl 6 S, 6 (j]ni) — 8 In 3 (in j) and the expression for the correlation length in the radiation dominated universe l|4.53|l . 

As in Minkowki space-time, the presence of a high conductivity plasma severely hinders the generation of magnetic 
fields. However, a noteworthy aspect is that up to the non-linear time the magnetic field is still correlated over the 
size of the scalar field domains rather than the diffusion length t;dif f ~ \A?/ctr. The diffusion scale determines the 
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spatial size of the region in which magnetic fields are correlated in the absence of non-equilibrium generation. The 
ratio between the domain size £ (77) given by 1)4. 53[) and the diffusion length scale £<ji / / iv) is given by 




2 



Where we have used the relations l|4.3[l , (|4.4|) , l|4.20|l and 14.4811 . Thus an important conclusion of this study is that 
the magnetic fields generated via spinodal decomposition are correlated over regions comparable to the size of scalar 
field domains which are much larger than the diffusion scale. 

The spectrum for the electric field can be obtained from that of the magnetic field by simply replacing k T>c — > 2?c- 
In the soft regime and for time scales ^ < t) C |f we have T>c — —k 2 /a R whereas kT> c ~ k/<jR. Therefore the 
electric field spectrum is suppressed by a factor k 2 /a R with respect to the magnetic field, namely 

ST(r),k) = ^rS°«(v,k). (5.88) 

°R, 

Thus, in a high temperature plasma with large conductivity the non-equilibrium processes favor the generation of 
magnetic photons instead of electric photons, and again equipartition is not fulfilled. 

The energy density on large scales > L again can be computed in closed form in the limits L 3> CiVnl) an d 
L <C We find in the first case from ea. (|5.8(l . 

2 13 ir-f Na 1 
Ap B (Vni,L) = j— 4 2 — - —j ; L> £,(ri n i) . (5.89) 



We find for the opposite case, 



3 x 2 10 7r 3 N a 

A ^ L ^ = X 2 ^a 2 R e°(v nl ) ; L K< aVnl) • (5 ' 90) 



The ratio of the magnetic energy density on scales larger than L at the spinodal time and the total radiation energy 
given by the Stefan-Boltzman law p 7 = 7r 2 T'* ; /15 results, 

Ap B ( Vn i,L) 3x2 13 tt^ Na 1 

r(r]nl,L) = = -2 j— , 14 g , ry , , 4 ; i > CWni) • 5.91) 

This result is the same as for the Minkowski space-time (see eq.(7.53) in ref.[23), except for a multiplicative factor 
/2 6 £ 6 (f?ni) — 8 In 3 (in -J) and the expression for the correlation length in the radiation dominated universe l|4.53ll . 

The factor (LTr) -4 is purely dimensional and is ultimately the detcrminining factor for the strength of the generated 
magnetic fields on a given scale. These combinations are invariant under cosmological expansion and arc determined 
by the ratio of the scales of interest today (galactic) to the thermal wavelength (today) of the cosmic microwave 
background radiation at the Wien peak. In particular LTr ~ 10 25 for L ~ 1 Mpc(today) [see ea. l|4.16|l ]. 

It is clear that the production during this regime is extremely small, due to the large values of (LT) A and of the 
ratio a 2 R J p? . In order to obtain an estimate for the amplitude of the seed magnetic field, we consider the following 
set of parameters: A = 10~ 2 , a = 10~ 2 , p = 10 14 GeV, Tr = 10 16 GeV (corresponding to a critical temperature 
T c = 10 15 GeV). We then obtain, 

r(L = IMpc) - 1CT 157 . (5.92) 

Therefore, the amplitude of the magnetic field generated during th e sp inodal stage is completely negligible. This 
result is similar to the result obtained in Minkowski space-time in ref. |27j and is expected on the basis of dimensional 
analysis. 



C. Magnetogenesis from the scaling regime 

In the scaling regime 77 >> r\ n \ the spectrum of the magnetic field is given by (|5.82|1 with the mode functions in the 
scaling regime given by (|4.59|) . 

The final expression for the leading contribution, given by (|5.82|) reveals a noteworthy aspect. As we have argued 
above, the modes k of astrophysical relevance today, were well outside the horizon during the radiation dominated 
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era between reheating and the QCD phase transition. The mode functions (|5.82|l attain the largest amplitude at long 
times for x = qrj < 2 — 3, thus momenta in the polarization loop that are within the horizon lead to generation of 
magnetic fields with long-wavelengths well outside the horizon. This we believe, is an important mechanism, loop 
corrections lead to a coupling between modes inside the horizon with those outside. Thus in this manner, causal 
fluctuations can actually lead to the generation of fields with wavelengths much larger than the horizon. 
Since krj -C 1 the power spectrum ()5.82(l takes the following form using the scaling mode functions 1)4. 59|) . 

Sb(v, k ) = (1 + 2 «o) 2 — \ A o\ 2 dq / J 2 (qvi) Vi *7i / hiom) V2 dr) 2 I(q, k, 771, r/ 2 ) , (5.93) 



where we set the exponentials equal to unity in eq. I|5.82[) since krj <C 1 and k <C (Jr and 

+ 1 l-X 2 

I(q,k,r]i,r]2) = J dx + fc2 _ 2kqx) 2 J2 ^ q2 + fc2 2kqx v ^ ' q2 + fc2 2kqx m ^ " 



Using the summation theorem[4f 



J 2 (Vq 2 + k 2 - 2kqx V ) = 4{<l2+ q f k2 2 2fcgX) ^(/ + 2) J l+2 (q V ) J l+2 (kn) C 2 (x) 



where the Cf(x) are Gegenbauer polynomials. For krj <C 1 the I = terms dominate and we can use the small 
argument behaviour of the Bessel functions J 2 (kn) — ^(kr)) 2 [l + <D(k 2 rj 2 )]. We finally obtain, 

/(<Z,Mi,>72) = A J 2(wi) Hqm) [l + 0(k 2 V 2 )] . (5.94) 
1 

Inserting ea. l|5.94|l into ea. l|5.93|l yields 

AT £.2 /"Oc 1 

S B (V, k) = (1 + 2n ) 2 J j- |^o| 2 / 4 K t J 2 2 (^) - -MoT) Mm)] ~ (V - W)} 2 [1 + ^(fc 2 ^)] , (5-95) 

where we used the formula 

" z Jf (/3z) dz = V - [j 2 {(3y) - Mpy) J 3 (py)] . (5.96) 
Since r/ 3> ?7_r we can neglect the terms with r)R and we find, for k <C f? -1 

rsN k 2 

S B (v,k)=V — — ^H 4 R r 1 7 [l + 0{k 2 f 1 2 )] ; © = 48.61.... (5.97) 
where we used ea. (l4.61|) and we computed numerically the integral 



dx 



x 



J 2 2 (x) - Ji(a;) J 3 (a:)r = 0.0005295 . . 



This integral is dominated by the region x > 1, namely, by modes that are inside the horizon. From the estimate 
1)4.24(1 . the corrections O (k 2 ^^^)^ are truly negligible between reheating and the QCD phase transition. 

The dependence on the conformal time ~ if is a direct consequence of the scaling form of the solution for the mode 
functions. The strong time dependence is a consequence of the causal relaxation of the Goldstone fields, a result of 
the phase ordering kinetics that entails that the size of the domains grow with the horizon. 

There are several important features of the above result which are noteworthy: 

• i: The exponential associated with the diffusion length cancels out, a reflection that the long time behavior of 
the integrals above are dominated by the upper limit. Hence the final result for the spectrum does not feature 
the exponential suppression with the diffusion length. 

• ii: The result for the spectrum only depends on the initial amplitud Aq and initial occupation number no in the 
combination |Ao| 2 (l + 2uq) which is constrained by the sum rule (|4.60|l . Hence the final spectrum is insensitive 
to the initial conditions on the mode functions or occupations, which in principle carry information of the early 
history beginning from the inflationary stage. This is a consequence of the scaling solution being a fixed point 
of the dynamics of the scalar field 26), l28l |29| . 
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• iii: A noteworthy result is that superhorizon magnetic fields are generated by the non-equilibrium dynamics 
of modes inside but near the Hubble radius. This is a consequence of the polarization loop, wherein the 
propagators correspond to momenta q and \q + k\. The momenta k corresponding the wavevector (scale) of the 
magnetic field is such that the wavelength is larger than the Hubble radius, but the momenta q corresponding 
to the charged scalar field fluctuations are inside the horizon. The correlation length of the charged scalar 
field is of the order of the Hubble radius. Thus acausal, superhorizon magnetic fields are generated by loop 
effects. Furthermore, the processes that yield the leading contribution are pair production, pair annihilation 
and low energy bremsstrahlung. These processes are not allowed in equilibrium because of energy momentum 
conservation, but are allowed in a rapidly expanding cosmology and strongly out of equilibrium as is the case 
under study. 

• iv: The out of equilibrium dynamics for very long superhorizon wavelengths cannot be captured with a kinetic 
description as advocated in ref. |14|. since kinetic descriptions rely on slow time variations of the distributions. 
Clearly such description will not be justified during the early and intermediate stages of the evolution after 
the phase transition because of the rapid (exponential) time evolution of the mode functions and therefore the 
distribution functions. The correct scaling dynamics cannot be captured within the kinetic description either 
since it describes the evolution of coherent domains correlated on scales comparable to the Hubble radius. Thus 
the formulation presented here based on the Schwinger-Dyson equations provides the necessary calculational 
scheme to describe reliably from first principles the generation of magnetic fields during and after non-equilibrium 
phase transitions. 

In order to reveal the enhancement during the scaling regime in a more transparent manner, it is convenient to use 
the relations (|4.3|l . (|4.13() . I|4.40[l and the explicit expression for the conductivity l|4.63[l in the form 

aR = c [ a ,N]*I*. , c [a,N] = — ^- ~ 0(1). (5.98) 
Then, the ratio r(L, 77) for L 3> ?/ is given by 

240tt V a 3 ( n \ 4 / M* \ 3 . 

r ^=Nc^N][LT R f {7m ) \m) ■ (5 - 99) 

where T> is given in eq. I|5.97|l . We note that in the final result l|5.99|l there is no dependence on the reheating 
temperature but only on the scale of symmetry breaking /j, the temperature at the time 77 and the scalar and gauge 
couplings. This is expected since the non-equilibrium processes begin in earnest after the phase transition, local 
thermal equilibrium prevailed between the time of reheating and the phase transition. 

The dependence on the scalar self coupling cx 1/A 2 is a hallmark of the non-perturbative nature of the growth of 
unstable modes and spinodal decomposition, it is ubiquitous in the non-equilibrium dynamics of phase transitions|26l 

Large scale magnetogenesis is more efficient for large symmetry breaking scale since the larger the symmetry 
breaking scale, the longer lasts the scaling stage. 

Consider for instance the case in which the symmetry breaking scale \i ~ 10 13 Gev and A ~ a ~ 10~ 2 , corresponding 
to a critical temperature of order of a GUT scale T c ~ 10 15 GeV and suppose that the scaling regime lasts until the 
electroweak phase transition scale, i.e. 77 is such that T(r]) = Tew ~ 10 2 GeV. Then the factor 

4 / „ r \ 3 



/ /x y( M *X 

\VXT EW J \Tew) 



10 



100 



compensates for the factor {LTr) 5 . Taking N and <?* of the order of 10 (these values are taken as representative and 
they can be changed simply in the final expressions) we can write the expression for the ratio as 



L y 5 (TewV 



Therefore 



,™,L).10-3^_^_J . (5.100) 

s rN [ 10~ 34 at the EW transition 
riT^)^) ~ I 1Q _ 14 attheQCDtransition . (5.101) 

Thus if the scaling regime lasts until a time between the EW and the QCD phase transitions the amplitude of the 
large scale magnetic fields is within the range necessary to be amplified by some dynamo models. The amplitude of 
the seed magnetic field is strongly dependent on the duration of the scaling regime. We have only focused on a scaling 
regime terminating either at the EW or QCD phase transition since there will surely be new phenomena associated 
with these that must be included in the dynamics of magnetogenesis. 
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D. Discussion 



Validity of the approximations: There are two main approximations that were used to obtain the results 
quoted above, i) the long-wavelength approximation kn <C 1 and ii) the weak coupling approximation. We now 
provide an estimate of the reliability of both these approximations to establish the limit of validity of our results. 



i): Long- wavelength approximation: In order to reach our final result for the rate r(L, 77) we have explicitly 
used a series of approximations which are valid for long wavelengths but whose validity must be checked before 
we reach any conclusion regarding the spectrum at small scales. In particular we must address the limits of 
applicability of the result (|5.99(l . This result has been obtained by integrating the magnetic spectrum on scales 
< k < kmax with k rn ax = 2tt / ' L m i n ; the formula for the magnetic spectrum was valid in the limit 

kmax Vmax ^ 1 • (5.102) 

In order to provide an estimate may take for rj max to be the (conformal) time at which the EW phase transition 
occurs, namely t\ew ~ 1 GeV~ . As discussed in the introduction, we are considering a situation in which the 
magnetic field is considered as a perturbation of a pre-existing thermal blackbody background. For consistency 
this requires that, 

r(L min ,r)wE) < 1 • (5.103) 

This relation translates in a condition 



kmax T]max < ^~~ 



CNa 3 ( T c V ( Tj 



- EW 



-i -1/5 



(5.104) 



^ 24 2 c 2 \Tew J V M * 

For T c ~ 10 15 GeV this gives k max r\ m ax <S 0.0176 which in turns translate into 

L > L min ~ 70 fm . (5.105) 

However, this is the comoving length normalized at the reheating time. In order to convert to the present time, 
we have to take in account the redshift 

Zji = ^~ 4 x 10 28 for T R ~ 10 15 Gev ; (5.106) 
To 

this gives 

Lmin\ toda y ~ 0.1 P C . (5.107) 

Thus, the approximations invoked are reliable to estimate the amplitude of primordial seeds on galactic scales 
or larger, today. 



ii) : Weak electromagnetic coupling: In order to study the amplitude for much smaller scales the calculations 
must be done without the long wavelength approximations invoked above. In this case we must expect a 
breakdown of perturbation theory and we cannot give a reliable estimate in the present framework. Furthermore, 
for scales well inside the Hubble radius, microphysical processes not included in our approximations, such as 
scattering between charged fields and between charged and gauge fields must be included. These processes 
will restore equilibrium between the different fields, if there is a substantial transfer of power from the charged 
fluctuations to the radiation field, this may lead to a change in the equation of state and the full backreaction 
on the metric must be included. At longer time scales the effects of the backreaction of the gauge fields on the 
dynamics of the scalar field, as well as the non-equilibrium contributions to equation of state and the Friedmann 
equations must be included self-consistently. 

Generation on short distance scales: For scales well inside the horizon during the scaling regime, namely 
qrj 1, we must account for causal microscopic processes that tend to equilibrate the electromagnetic fields 
generated by the non-equilibrium processes. In order to understand these processes we must look at the kinetics 
of equilibration. The mode functions for wavevectors well inside the horizon are Minkowski-like, of the form 
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where the coefficients a q , f3 q must be determined from a full numerical evolution. However the constancy of the 
Wronskian entails that 

KlM/3 9 l 2 = i 

which suggests the identification \a q \ = 1 + Af q ; \fl q \ = Af q , N q is the number of (asymptotic) quanta created 
during the time evolution. This form of the asymptotic mode functions leads to the equipartition between the 
electric and magnetic field generation, since spatial and time derivatives are the same. In turn this entails that 
we can understand the generation of electric and magnetic fields by obtaining a kinetic equation for the number 
of photons. Such kinetic equation must necessarily be of the form 

= [i + N k (r,)]T>(r,) - N k ( V )T<( v ) 

which displays the familiar gain minus loss contributions in terms of the forward and inverse rates. Eventually 
a steady state will be reached which will describe a stationary distribution of photons. The computation of 
the forward [r^(ry)] and inverse [rj5(?y)] require a detailed knowledge of the distribution A/",,|46j| since these 
generalized rates are functionals of these occupation numbers. Clearly such computation lies beyond the scope 
of this article and is a task that we will undertake elsewhere. However, the kinetic equations above will tend to 
an equilibrated state of local thermodynamic equilibrium. 

• Effect on the LSS: It is important to estimate the effect of the magnetic field on scales corresponding to 
those of the last scattering surface, which today are L^ss ~ 100 Mpc. From eq. I|5.101|l we see that at the 
electroweak temperature t{Tew ■> Llss) ~ 10 -44 , taking the fourth root we can provide an estimate of the 
temperature fluctuation induced by the primordial magnetic field ST/T\ LSS ~ 

[r{T EW ,L LSS )]* ~ 10- 11 which 

is negligible compared to the CMB temperature fluctuation at this scale ~ 10 -5 . On the other hand, a similar 
estimate at the time of the QCD phase transition gives 8T/T\ LSS ~ 10~ 6 which is marginally compatible with 
the current observations. Thus the reliability of the approximation of weak gauge coupling combined with the 
effects on the temperature anisotropy at the last scattering surface seem to lead us to conclude that if a phase 
transition during a radiation dominated era occurs near the GUT scale and results in a scaling stage, our results 
for primordial magnetogenesis will be reliable down to the scale of electroweak symmetry breaking. 



VI. CONCLUSIONS 



Wc review here large scale primordial magnetogenesis during a phase transition in the radiation dominated era after 
reheating in a model of TV-charged scalars coupled to an abelian gauge field |27|. The spectrum of the magnetic field 
generated during the non-equilibrium evolution was computed using the formulation recently introduced in ref.|27|. 
The dissipative effects of the conductivity are included by separating the contribution from hard modes (with momenta 
of order T) to the polarization tensor of the gauge fields. These modes are always in local thermodynamic equilibrium. 
The non-perturbative, non-equilibrium dynamics of the scalar field after the phase transition was studied in the large 
N limit. The dynamics after the phase transition features two distinct stages: an early and intermediate time, spinodal 
stage, which is dominated by the growth of long-wavelength fluctuations, followed by a scaling regime during which the 
scalar field becomes correlated over horizon-sized domains. During both regimes, strong non-equilibrium fluctuations 
lead to large current-current correlation functions which entail the generation of magnetic fields. The scaling regime 
is the most effective for primordial magnetogenesis since this stage lasts the longest. During this stage magnetic 
fields with superhorizon wavelengths are generated via loop effects, the scalar field momenta in the polarization loop 
corresponds to wavelengths of the order of or shorter than the horizon. Thus causal scalar field fluctuations lead to the 
generation of magnetic fields on superhorizon scales. In particular the processes that yield the leading contribution 
to magnetic field generation are i): pair production, ii) pair annihilation and iii) low energy bremsstrahlung, which 
cannot occur in equilibrium because of energy momentum conservation, but do occur out of equilibrium because of 
the rapid time evolution due to the expansion and the evolution of the scalar field out of equilibrium. 

The generation of magnetic field is hindered by the large conductivity of the plasma and equipartition between 
electric and magnetic fields does not hold. The spectrum of the primordial magnetic field is insensitive to the 
magnetic diffusion length which is sub-horizon during the radiation era. 

Our final result for the spectrum generated during the scaling regime is given by eq. I|5.97|l . The ratio of the energy 
density of the magnetic fields on scales larger than L to the energy density in the cosmic background radiation 
r(L,rj) = pb(L,t}) / p cm b{L,rf) is given by equation (|5.99() . For values of N, and the gauge coupling consistent with 
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particle physics models we find that 
Therefore, 

> n J 10 -34 at the EW transition , . 

r{I {r)),L) ~ | 1Q _ 14 attheQCD transition . ^) 

Therefore, the large scale primordial magnetic fields generated during the scaling stage after a phase transition may 
be a plausible mechanism to generate primordial magnetic fields which will be further amplified by the collapse of 
protogalaxies and by astrophysical dynamos. 

Probably a phase transition at a temperature much larger than the electroweak leading to a scaling regime lasting 
until the QCD phase transition is ruled out by the temperature inhomogeneities at the last scattering surface. Further- 
more the generation of electromagnetic fields on sub-horizon scales requires a full kinetic equation that incorporates 
the microscopic causal processes that lead to thermalization, the study of these is beyond the scope of this article. 

Magnetogenesis after the QCD phase transition? The model that we studied here is assumed to describe 
the robust features from the non-equilibrium dynamics of a charged sector coupled to a (hyper) charge gauge field. 
GUT's or SUSY theories may provide the corresponding framework. 

However we now argue that precisely the model studied here can actually describe the non-equilibrium dynamics 
after the QCD phase transition(s). After hadronization and chiral symmetry breaking most of the hadrons produced 
will be pions, at least this is the experimental situation in ultrarelativistic heavy ion collisions. Neglecting the charge 
form factor (which is justified for momenta much smaller than the p meson mass m p ~ 770 Mev) the charged pions 
couple to the electromagnetic field with minimal coupling. The chiral transition is conjectured to be in the same 
universality class as the 0(4) linear sigma model 47]. Thus the model presented in this article, is the low energy 
effective field theory for the triplet of pions, two charged and one neutral. 

Thus we conjecture that the study in this article can be applied to the study of the generation of magnetic 
fields by long-wavelength pions. Therefore the analysis of this article can apply to magnetogenesis during the chiral 
phase transition in QCD. While the charged pions couple to electromagnetism via the minimal coupling in the long- 
wavelength limit, the neutral pion couples to the electromagnetic field through the chiral anomaly 7r° — > 2j and such 
process will also produce magnetic fields. 

Since the ratio of the relaxation rate via the strong interactions to the expansion rate 



TqcD r ~ 10 28 (6.3) 

H(T]QCD) 

then if the chiral phase transition is second order, it will occur in LTE. However if it is first order as suggested by 
the fact that the pions are massive, then non-equilibrium effects may arise and primordial magnetic fields may be 
generated by charged pion bremsstrahlung or pair annihilation or by neutral pion decay. These possibilities can be 
studied with the linear sigma model effective field theory. 

Acknowledgements: We would like to thank Norma Sanchez for organizing this wonderful and very stimulating 
IX Chalonge School of Astrofundamcntal Physics. The authors thank M. Giovannini, for useful discussions. D. B. 
and M. S. thank N.S.F. for support through grants PHY-9988720 and NSF-INT-9815064. 



[1] E. N. Parker, Cosmical Magnetic fields, Clarendon, Oxford, 1979. 

Ya. B. Zel'dovich, A. A. Ruzmikin, D.D. Sokoloff, Magnetic Fields in Astrophysics, Gordon and Breach, New York, 1983. 
[2] P. P. Kronberg, Rept. Prog. Phys. 57, 325 (1994). 
[3] D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 (2001). 
[4] L. Widrow, Revs, of Mod. Phys. 74, 775 (2002). 

[5] M. Giovannini, hep-ph/0208152 to appear in the Proceedings of the 7th. Colloque Cosmologie, H J de Vega and N. Sanchez, 

editors. M. Giovannini, hep-ph/0111220 

[6] For a recent review on observations see: J.-L. Han and R. Wielebinski, astro-ph/0209090 

[7] D. Grasso, hep-ph/0002197 P. Elmfors, K. Enqvist and K. Kainulainen, Phys. Lett. B440, 269 (1998); M. Giovannini and 

M. E. Shaposhnikov, Phys. Rev. D57, 2186 (1988). 
[8] T. Kahniashvili, A. Kosowsky, A. Mack and R. Durrer, astro-ph/0011095 (2000); A. Mack, T. Kahniashvili, A. Kosowsky, 

Phys. Rev. D65 123004 (2002); R. Durrer, P.G Ferreira, T. Kahniashvili, Phys. Rev. D61, 043001 (2000); J. Adams, U. 

H. Danielsson, D. Grasso, H. Rubinstein, Phys. Lett. B388, 253 (1996). 



35 



[9] A. D. Dolgov, hep-ph/0110293 (2001). 
[10] C. J. Hogan, Phys. Rev. Lett. 51, 1488 (1983). 
[11] M. S. Turner and L. M. Widrow, Phys. Rev. D37, 2743 (1988). 

[12] W. D. Garretson, G. B. Field, S. M. Carroll, Phys. Rev. D46, 5346 (1992 ); B. Ratra, Ap. J . 391, LI, (1992). 
[13] M. Giovannini and M. E. Shaposhnikov Phys. Rev. D 62 , 103512 (2000) hep-ph/0004269 . 
[14] M. Giovannini and M. E. Shaposhnikov, hep-ph/0011105 (2000) 

[15] A. C. Davis, K. Dimopoulos, T. Prokopec and O. Tornkvist, Phys. Lett. B 501, 165 (2001); K. Dimopoulos, T. Prokopec, 
O. Tornkvist and A. C. Davis, Phys. Rev. D 65, 063505 (2002); T. Proko pec, |astro -ph/0 106247] 

[16] E. A. Calzetta and A. Kandus, Phys. Rev. D 65, 063004 (2002) astro-ph/0110341 . E. A. Calzetta, A. Kandus and 
F. D. Mazzitelli, Phys. Rev. D 57, 7139 (1998) astro-ph/9707220 . 

[17] B. A. Bassett, G. Pollifrone, S. Tsujikawa, F. Viniegra; Phys. Rev. D63 103515 (2001). 

[18] F. Finelli and A. Gruppuso, Phys. Lett.B502 216 (2001); F. Finelli, astro-ph/0007290 

[19] D. S. Lee, W. 1. Lee and K. W. Ng, Phys. Lett. B542, 1 (2002). 

[20] A. Dolgov, Phys. Rev. D 48, 2499 (1993). 

[21] A.-C. Davis and K. Dimopoulos, Phys. Rev. D55 7398, (1997). 
[22] T. Vachaspati, Phys. Lett. B 265, 258 (1991). 

[23] K. Enqvist and P. Olesen, Phys. Lett. B 319, 178 (1993); ibid 329, 195 (1994); K. Enqvist, in Strong and Electroweak 

Matter '97, Ed. F. Csikor, Z. Fodor, World Scientific, Singapore, 1998. 
[24] G. Baym, D. Bodeker, L. McLerran, Phys. Rev. D53 662 (1996). 
[25] J. Ahonen and K. Enqvist, Phys. Rev. D57, 664 (1998). 

[26] D. Boyanovsky and H. J. de Vega, Phys. Rev. D 61, 105014 (2000); D. Boyanovsky, H. J. de Vega and R. Holman, 
Non-equilibrium phase transitions in condensed matter and cosmology: spinodal decomposition, condensates and defects, 
in: Topological Defects and the Non-Equilibrium Dynamics of Symmetry Breaking Phase Transitions, Eds. Y. M. Bunkov 
and H. Godlrin, Nato Science Series, Kluwer, 2000. 

[27] D. Boyanovsky, H.J. de Vega and M. Simionato, Magnetic field generation from non- equilibrium phase transitions, 
hep-ph/0208272 to appear in Phys. Rev. D; Large scale magnetogenesis from a non- equilibrium phase transition in the 
radiation dominated era, hep-ph/0211022 

[28] N. Turok and D. N. Spergel, Phys. Rev. Lett. 66, 3093 (1991); D. N. Spergel, N. Turok, W. H. Press and B. S. Ryden, 
Phys. Rev. D43, 1038 (1991). 

[29] R. Durrer, M. Kunz, A. Melchiorri, Phys. Rep. 364 1 (2002); R. Durrer, M. Kunz, A. Melchiorri, Phys. Rev. D59 123005 

(1999) . 

[30] D. Boyanovsky, D.-S. Lee and A. Singh, Phys. Rev. D48, 800, (1993); see also ref. |H. 

[31] D. Boyanovsky, H. J. de Vega and M. Simionato, Phys.Rev. D63, 045007 (2001); D. Boyanovsky and H. J. de Vega, 
Phys.Rev. D65, 085038 (2002). 

[32] D. Boyanovsky, D.-S. Lee and A. Singh, Phys. Rev. D48, 800, (1993); D. Boyanovsky, H. J. de Vega, R. Holman, 
D.S-Lee and A. Singh, Phys. Rev. D51, 4419 (1995), D. Boyanovsky, D. Cormier, H.J. de Vega and R. Holman, Phys. 
Rev D55, 3373 (1997); D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman and S. Prem Kumar, Phys. Rev. D57, 
2166 (1998); D. Boyanovsky, H. J. de Vega and R. Holman, Proceedings of the Second Paris Cosmology Colloquium, 
Observatoire de Paris, June 1994, pp. 127-215, H. J. de Vega and N. Sanchez, Editors, World Scientific, 1995; Advances 
in Astrofundamental Physics, Erice Chalonge Course, N. Sanchez and A. Zichichi Editors, World Scientific, 1995; D. 
Boyanovsky, H. J. de Vega and R. Holman, p. 183-270 in the Proceedings of the Vth. Erice School 'D. Chalonge', 'Current 
Topics in Astrofundamental Physics', N. Sanchez and A. Zichichi, editors, World Scientific. D. Boyanovsky, D. Cormier, 
H. J. de Vega, R. Holman and S. Prem Kumar, Erice Chalonge School, Current Topics in Astrofundamental Physics, 
N. Sanchez and A. Zichichi Editors, p. 119-194, Kluwer, 1998; D. Boyanovsky and H. J. de Vega, Out of equilibrium 
fields in self consistent inflationary dynamics: Density fluctuations astro- ph/0006446 in the Proceedings of the Vllth. Erice 
Chalonge School on Astrofundamental Physics, p. 37, N. Sanchez ed., Kluwer, Series C, 2000. 

[33] D. Boyanovsky, H. J. de Vega, R. Holman, J. F. J. Salgado, Phys. Rev. D54, 7570 (1996); D. Boyanovsky, C. Destri, H. 
J. de Vega, R. Holman and J. F. J. Salgado, Phys. Rev. D57, 7388 (1998); D. Boyanovsky, H. J. de Vega, R. Holman and 
J. F. J. Salgado, Phys. Rev. D59, 125009 (1999). 

[34] D. Boyanovsky, D. Brahm, R. Holman and D. S. Lee, Phys. Rev. D54, 1763 (1996); D. Boyanovsky, H. de Vega and S.-Y. 
Wang, Phys. Rev. D61, 065006 (2000). 

[35] D. Boyanovsky, H. J. de Vega, R. Holman and M. Simionato, Phys. Rev. D60 (1999) 065003; D. Boyanovsky and H. J. de 
Vega, Phys. Rev. D59 (1999) 105019; S.-Y. Wang, D. Boyanovsky, H.J. de Vega, and D.-S. Lee, Phys. Rev. D 62, 105026 

(2000) . 

[36] G. Baym and H. Heiselberg, Phys. Rev. D56, 5254 (1997). 
[37] P. Arnold, G. D. Moore, L. G. Yaffe, JHEP 0011, 1 (2000). 

[38] D. Boyanovsky, H.J. de Vega and R. Holman, Phys. Rev. D49, 2769 (1994); D. Boyanovsky, D. Cormier, H. J. de Vega, 

R. Holman, A. Singh, M. Srednicki, Phys. Rev. D56,1939 (1997). 
[39] J. Schwinger, J. Math. Phys. 2, 407 (1961); K. T. Mahanthappa, Phys. Rev. 126, 329 (1962); P. M. Bakshi and K. T. 

Mahanthappa, J. Math. Phys. 41, 12 (1963); L. V. Keldysh, JETP 20, 1018 (1965); K. Chou, Z. Su, B. Hao And L. Yu, 

Phys. Rep. 118, 1 (1985); A. Niemi and G. Semenoff, Ann. of Phys. (NY) 152, 105 (1984); N. P. Landsmann and C. G. 

van Weert, Phys. Rep. 145, 141 (1987); E. Calzetta and B. L. Hu, Phys. Rev. D41, 495 (1990); ibid D37, 2838 (1990); J. 

P. Paz, Phys. Rev. D41, 1054 (1990); ibid D42, 529(1990). 
[40] D. Boyanovsky, H. J. de Vega and R. Holman, Proceedings of the Second Paris Cosmology Colloquium, Observatoire de 



36 



Paris, June 1994, pp. 127-215, H. J. de Vega and N. Sanchez, Editors, World Scientific, 1995; Advances in Astrofundamental 
Physics, Erice Chalonge Course, N. Sanchez and A. Zichichi Editors, World Scientific, 1995. 

[41] See for example: D. Boyanovsky, H. J. de Vega and R. Holman, p. 183-270 in the Proceedings of the Vth. Erice School 
'D. Chalonge', 'Current Topics in Astrofundamental Physics', N. Sanchez and A. Zichichi, editors, World Scientific. D. 
Boyanovsky, D. Cormier, H. J. de Vega, R. Holman and S. Prem Kumar, Erice Chalonge School, Current Topics in 
Astrofundamental Physics, N. Sanchez and A. Zichichi Editors, p. 119-194, Kluwer, 1998 and references therein. 

[42] J. I. Kapusta, Finite Temperature Field Theory, Cambridge Monographs on Mathematical Physics, Cambridge Univ. Press, 
1989. 

[43] M. Le Bellac, Thermal Field Theory, Cambridge University Press, 1996. 
[44] A. Dolgov and J. Silk, Phys. Rev. D 47, 3144 (1993). 

[45] I. S. Gradshteyn and I. M. Ryshik, Table of Integrals, Series and Products, Academic Press, 1980; N. N. Lebedev, Special 

Functions and their Applications, Prentice Hall, 1965. 
[46] D. Boyanovsky, H. J. de Vega and M. Simionato, Phys. Rev. D 61, 085007 (2000). 

[47] K. Rajagopal and F. Wilczek, Nucl. Phys. B399, 395 (1993); Nucl. Phys. B404, 577 (1993); K. Rajagopal, in 'Quark-Gluon 
Plasma 2', Ed. by R. C. Hwa (World Scientific, Singapore) (1995). 



